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ABSTRACT 

A review  of  the  work  on  distributions  of  characteristic  (ch.J  roots  in 
real  Gaussian  multivariate  analysis  has  been  attempted  surveying  the  develop- 
ments in  the  field  from  the  start  covering  about  fifty  years.  The  exact  null 
and  non-null  distributions  of  the  eh.  roots  have  been  reviewed  and  sub-asymptotic 
and  asymptotic  expansions  of  the  distributions  mostly  for  large  sample  sites 
studied  by  various  authors,  have  been  briefly  discussed.  Such  distributional 
ies  of  four  test  criteria  and  a few  loss  important  ones  which  are  functions 
h.  ***>ots  have  further  been  discussed  in  view  of  the  power  comparisons  made 
in  coi  ion  wich  tests  of  three  multivariate  hypotheses.  In  addition,  one- 
sample  ease  has  also  been  considered  in  terms  of  distributional  aspects  of  the 
eh.  roots  and  criteria  for  tests  of  two  hypotheses  on  the  covariance  m w.-ix. 

A brief  critical  review  has  also  been  attempted.  For  convenience  in  organiza- 
tion, the  review  has  been  given  in  two  parts:  Part  I.  Null  distributions  and 

Part  II.  Non-null  distributions. 


Invited  paper  _ The  Canadian  Journal  of  Statistics 

’This  work  was  supported  in  part  by  the  Office  of  Naval  Research  contract 
N00014-75-C-0455  at  Purdue  University.  Reproduction  in  whole  or  part  is 
permitted  for  any  purpose  of  United  States  Government. 


PART  I 

NULL  DISTRIBUTIONS 


1.  INTRODUCTION 

In  tlu*<  part  of  the  review,  the  null  distribution  of  the  eh.  roots  of 
sample  matrices  arising  in  thuo  tests  of  hypotheses  will  he  discussed  first 
and  then  the  one-snnple  case  (Section  2J.  In  Section  3,  four  criteria,  namely, 
l!  Hotel  line's  trace,  2)  PillaiS  trace,  3)  Key's  largest  (smallest)  root  and 
U Wilks'  criterion,  recommended  for  tests  of  the  three  hypotheses  and  their 
optimal  properties,  and  further  three  other  less  interesting  criteria  will  be 
briefly  treated.  In  addition,  tests  of  two  hypotheses  on  the  covariance  matrix 
will  he  considered  in  the  one-sample  case  and  also  some  tost  statistics  in  this 
connection.  Section  *1  consists  of  the  study  of  the  densities  of  individual 
ch.  roots  for  which  four  different  methods  are  treated.  Distribution  problems 
of  the  two  trace  criteria  are  studied  in  Section  S discussing  exact,  approxi- 
mate and  asymptotic  cases  and  showing  the  relation  between  the  moments  in  the 
two  cases.  Section  6 contains  a treatment  of  the  distribution  of  Wilks' 
criterion,  the  multivariate  beta  distribution  and  associated  independent  beta 
variables.  The  exact  and  approximate  distributions  are  considered.  Other 
statistics  of  less  importance  are  discussed  in  Section  T^namely,  Wilks'  statis- 
tic, Pi llni's  harmonic  mean  criteria  and  Bagai's  statistic.  Further,  the  one- 
sample  case  is  taken  up,  treating  the  distribution  problem  of  the  likelihood 
ratio  statistic  for  testing  the  covariance  matrix  equal  to  a specific  matrix 
and  those  of  sphericity,  Wilks'  sample  generalized  variance,  ratios  of  ch. 
roots  and  elementary  symmetric  functions  of  the  roots. 

The  review  does  not  cover  the  material  on  complex  Gaussian  distributions 
although  they  arc  treuted  in  a large  number  of  papers  listed  in  the  References 
(For  example,  see  James,  1964,  Krislmaiah,  1976). 
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2.  NULL  DISTRIBUTION  OP  THL  CHARACTERISTIC  ROOTS 

In  multivariate  statistical  analysis  often  we  wish  to  test  one  or  more  of 
the  following  three  hypotheses:  (Pillnl,  1957,  1960) 

1 1 Equality  of  covariance  matrices  of  two  p-varlate  normal  populations; 

II)  Equality  of  the  p-dinensional  mean  vectors  of  % p-varlate  normal 
populations  having  a common  covariance  matrix  known  ns  MANOVA  (alternately 
general  linear  hypothesis) ; and 

III)  Independence  between  a p-set  and  a q-set  (p  ^ q)  In  a (p+q)-variate 
noir.nl  population  Known  as  zero  canonical  corrclntion  (2CC). 

Tests  proposed  for  the  above  hypotheses  are  generally  invariant  tests 
(Roy,  195?,  Lehmann,  1959)  which  have  been  shown  to  depend,  under  the  null 
hypotheses,  only  on  the  characteristic  roots  of  matrices  based  on  sample 
observations.  lor  example,  in  1),  the  tests  depend  on  the  characteristic 
roots  of  s j fSj  + S,)"1,  where  S^  and  S,  denote  the  sum  of  product  (SP)  matrices 

with  n,  and  n,  degrees  of  freedom  (df)  and  where  both  ?.re  almost  everywhere 

* -1 
positive  definite  (aepd) . Tlius  S,(^^  + S2)  is  acpd  whence  it  follows  that 

all  the  p ch.  roots  arc  greater  than  zero  and  less  than  unity,  in  II),  the 
♦ * ' + 

matrix  is  S (S  + S)  , where  .S  denotes  the  between  SP  matrix  of  means 
«•  «.  ■*  M 

weighted  by  the  sample  sizes  with  £-1  df  and  S denotes  the  within  SP  matrix 

■* 

(pooled  from  the  SP  matrices  of  * samples)  with  N-fc  df  where  N is  the  total 

of  the  sices  of  the  i samples  where  S is  acpd  and  S is  at  least  positive 

**  (* 

semi -definite  of  rank  i = min(p,fc-l).  Thus,  ac,  s of  the  ch.  roots  arc 
greater  than  sera  and  less  than  unity  and  p-s  remaining  roots  arc  aero. 

In  III),  the  matrix  is  S^jSpS^Sp,  where  Sp  is  the  SP  matrix  of  the  sample 
of  observations  on  the  p-set  of  variables,  S2,,  that  on  the  q-set  and  Sp, 
the  SP  matrix  between  the  observations  on  the  p-set  and  those  on  the  q-set. 

If  p P+R  < n'»  the  sample  size,  then  ac  the  p ch.  roots  of  this  matrix 
arc  greater  than  zero  and  loss  than  unity. 

In  each  of  the  three  cases  above,  if  the  hypotheses  to  be  tested  is  true, 
the  s < p non-zero  ch.  roots  b.,i=l,...,s,  where  0 < b,  < b,  <..<  b <1, 
have  the  same  form  of  the  joint  density  function  which  was  independently 
obtained  in  1939  by  the  five  authors,  Fisher  (1939),  Girshick  (1939),  Hsu 
(1939),  Mood  (1951)  and  Roy  (1939).  This  joint  density  can  be  written  in  the 
form 
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’•t  i{n-i,)n  .1  C>  -h.»,  n < h.  ...< 
.-i  1 4 » -i  11  1 


l»  - 1, 


where 

C-.2I  C(s,m,nl  * si  (r[${2iatdn»*M«dl|/lM$(2ttM*n)r(f(2nM*!))r(l/2)J. 

Here  m and  n are  to  be  understood  differently  for  different  situations.  For 
example,  in  I), 


(2.3)  m " f (nj-p-1),  n ■ f.(n,-p-l). 

In  II), 

(2.4)  » - MiU-U-pl-1),  n - M(N-t)-p-l]* 

In  III, 

12.5)  m ■ f.(q-p-l) , n ■ ^[(n’-l-q)-p-l). 

Note  that  n and  n in  (2.3)  to  (2.5)  could  written  in  terns  of  Vj  and  v,  df, 
for  example,  in  I),  Vj»nj,v,«n,;  in  II),  Vj«l-1  and  in  III),  Vj»q  and 

-1  1 

v,»n'-l-q,  where  q is  the  of  of  the  matrix  S15S,,S,,  and  n'-l-q,  of 
-1  * 

sn"VS”SP*  Further,  (2.2)  may  be  rewritten  in  terms  of  Vj  and  v,. 


Alternately,  there  is  interest  in  considering  the  eh.  roots  of  S.S'1, 

S S or  (^12^22^12^ ^ll"!!l2^22^1  2^  given  by  * b^/O-b^),  i*l,...,$, 
where  0 < f,  fc  < ».  The  joint  density  of  fj,...,fs  may  be  obtained 


1 " s 
from  (2.1)  in  the  form 


m*n*s*l 

(2.6)  g(f,,...,f)  - c(s,»,n)  n (fT/O+fJ)  n (f,-f.),  o < f.  *..<f( 

i«l  i>j  J 


< •. 


The  form  of  the  density  in  (2.1)  will  be  called  Type  I and  that  in  (2.6), 
Type  II. 

Further,  if  « nb.,  i»l,..,s,  as  n-^»,  the  density  of  gj,...,gs  is 
given  b'*,  (Nanda,  1943b)  s 

(2.7)  h(g.,...,g  ) « K(s,m)  n (g?  e x)  n (g.-g.),  0 < g.  <..<  g < •, 
s i.i  1 i>)  x 3 l - - s 


as  s 

where  K(s,m)  ■»  nJ;5/  n (r[f.(2ra*i+l)]r(i/2) }. 

i-1 

The  density  (2.7)  is  of  interest,  especially  under  the  hypothesis  I) 
since  has  a Mishart  distribution  (Wishart,  1928)  with  df  and 
covariance  matrix  Ej,  denoted  by  KCp.Vj.Ej),  the  density  of  ch.  roots  of 

SjEj^  when  Ej  is  '..nown,  could  be  written  in  the  above  form. 


Again,  it  is  easy  to  sec  that  many  of  the  commonly  used  distributions  in 
statistics  are  special  cases  of  (2.1) for  (2.6)  or  (2.7)).  For  example,  if  p«l 
in  (2.3),  F distribution  is  obtained  from  (2.1)  or  more  directly  from  (2.6)  and 
chi -squared  distribution  from  (2.7);  if  *»2  in  '2.4),  similarly  is  obtained 
from  (2.6)  the  density  of  Mahalanobis1  D“  (Mahalanobis,  1930,  1936), 

(Hotelling's  T“  (Hotelling,  1931).),  if  A»2  and  p-1  Fisher's  z (student's  t); 
if  pi  in  (2.5 j , from  (2.1)  is  obtained  the  density  of  the  multiple  correlation 
coefficient,  if  p»l,  q*l  that  of  the  simple  correlation  coefficient;  and  others. 

An  outline  of  the  derivation  of  the  distributions  given  above  will  be 
deferred  to  later  sections  (Section  6;  Section  9,  Lemma  9.1).  For  the 
derivation,  in  addition  to  the  r*.  ‘crenecs  preceding  (2.1),  references  may  be 
made  to  James  (1954),  Dccmer  and  Olkln  (1951)  and  Olkin  and  Roy  (1954).  In 
order  to  discuss  the  usefulness  of  these  distributions  it  may  be  appropriate 
to  consider  briefly  the  important  test  statistics  in  the  literature  which  are 
functions  of  the  ch.  roots.  These  will  be  considered  in  the  next  section. 

3.  SOME  TOST  CRITERIA 

Four  statistics  which  are  commonly  used  for  tests  of  the  three  hypotheses 
considered  in  the  previous  section  are  as  follows: 

1)  Hotelling's  trace,  ■ £ f . , (Pillai,  1954,  1955),  with  Hotelling's 

i-1  1 

To  * v->ufS)»  (Hotelling,  1344,  1947,  1951); 

2)  Pillai 's  trace,  V(s)  » l b.,  (1954,  1955); 

i»l  1 

3)  Roy's  largest  (smallest)  root,  (1945,  1957);  and 

4)  Wilks'  criterion,  W(s)  ■ R (1-b.),  (1932). 

i*l  1 

The  statistic  1)  is  known  as  Hotelling's  trace  and  less  frequently  as 
Lawley-Hotelling  trace.  Tq  which  was  proposed  by  Hotelling  in  1944  as  a 
generalized  T-tost  and  measure  of, multivariate  dispersion  has  been  previously 
considered  by  several  authors,  lawley  (1938)  as  a generalization  of  Fisher's 
z-test,  Bartlett  (1939),  and  Hsu  (1940a) who  also  obtained  the  null  density  for 
s*2  and  the  first  and  second  non-central  moments  of  the  statistic  explicitly. 
The  statistic  2)  known  as  Pillai's  trace  which  was  proposed  by  Pillai  in  1954, 
(Pillai,  1954,  1955)  was  considered  by  Bartlett  (1939),  Hotelling  (1947)  and 
Nanda  (1950).  Wilks'  criterion  proposed  in  1932  for  test  of  hypothesis  II) 
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is  some tines  ttfiiuM’il  by  Kjlk.'  A ( older  mu.it  ion),  IJ  , „ , and  more 
tsl  P ,J1,V2 

recent  l>  by  W (see  I’ll  In  i ami  .luyuchnmlr.in  . 1 *»f*7 , 1968)  and  is  the  oldest 

of  the  four  test  criteria  hut  although  the  exact  distribution  was  studied  and 

percentage  points  obtained  by  Schntzoff  (1966b) , Pi  1 lai  and  Gupta  (1969)  and 

Lee  (1972)  only  three  decades  lnter,a  chi-squared  approximation  to  its  density 

suggested  by  Bartlett  (1938,  1947a)  and  a beta  approximation  by  Rao  (1951, 

1952)  enabled  the  early  use  of  the  criterion,  further,  Pillai  (1955)  proposed 

Kilkf.'  criterion  for  test  of  hypothesis  1)  as  well.  The  use  of  the  other 

three  criteria  was  facilitated  by  the  pioneering  statistical  tables  of  Pillai 

(1957,  I960).  The  monotonicity  of  powers*  of  these  criteria  with  respect  to 

each  population  deviation  parameter  (see  Section  9)  has  been  studied  by 

various  authors  (Das  Gupta,  Anderson  and  Mudholknr,  1964  , Anderson  and 

Das  Gupta,  1964a,  1964b  , Roy  and  Mikhail,  1961  , Mikhail,  1962  , Srivastav?, 

1964  , Mudholkar,  1965  , Perlman,  1974  , Eaton  and  Perlman,  1974  ).  for  II), 

the  admissibility  of  l/5^  and  the  largest  root  has  been  established  by  Ghosh 

(1964)  against  unrestricted  alternatives  and  Schwartz  (1964)  that  of  in 

the  same  sense.  Kiefer  and  Schwartz  (1965)  have  shown  that  test  is 

admissible  Bayes,  fully  invariant,  similar  and  unbiased.  They  have  also 

shown  that  is  admissible  Bayes,  under  a restriction,  although 

admissibility  could  be  established  without  this  restriction.  Giri  (1968) 

has  shown  that  is  locally  best  invariant  test  for  I)  against  one-sided 

alternatives.  Pillai  and  Jayachandran  (1967,  1968)  have  made  exact  power 

comparisons  in  the  two-roots  case  of  tests  of  each  of  the  three  hypotheses 

based  on  the  four  criteria  and  have  shown  that  and  compare 

favorably  and  behave  somewhat  in  the  same  manner  in  regard  to  the  three 

hypotheses  but  the  largest  root  has  lower  power  than  the  other  three  when 

there  are  more  non-zero  deviation  parameters  than  one.  Schatzoff  (1966a) 

has  similar  findings  in  the  case  of  II)  in  his  monte  carlo  study.  Again 

for  II),  Fujikoshi  (1970)  has  obtained  some  approximate  powers  for  l)^,V^  and 

and  Lee  (1971b)  some  approximate  powers  for  s*3  and  4,  and  Pillai  and 

Sudjana  (1974)  exact  powers  for  for  all  the  three  hypotheses  and 

approximate  powers  for  s=3,4  and  5 for  II). 


Although  such  aspects  should  be  discussed  under  non-central  distributions, 
this  discussion  is  included  here  to  give  some  idea  of  the  relative  importance 
of  the  different  test  criteria  which  would  help  the  understanding  of 
subsequent  sections. 


l.xaet  power  comparison  of  tests  of  I)  based  on  individual  roots  have  been 
made  bv  1‘ill.u  .uni  Al-Ani  f 1*170)  and  for  11}  and  III}  by  PilUl  and  Dotson  (1969) 
Although  the  largest  root  has  generally  more  power  than  the  other  roots,  in  small 
samples  it  appears  that  the  smallest  root  or  the  middle  root  can  have  more 
power  than  the  lurgcst. 

It  nay  be  noted  that  Wilks'  criterion  isafunction  of  the  likelihood  ratio 
criterion  (ire}  for  II)  and  111}  but  not  for  II.  For  I),  the  statistic 

n b.  (1-b.)  * is  a function  of  the  ire  modified  by  Bartlett  (1937) replacing 

i-1  1 1 

sample  site  by  df.  The  unbiasedness  of  the  test  has  been  shown  by  Suguira  and 
Nagao  (1908)  against  one-sided  or  two-sided  alternatives  while  Wilks'  criterion 
has  been  shown  so  far  to  have  optimum  properties  only  for  one-sided  alternatives. 
Again  for  two-sided  alternatives  Roy  (1957)  has  shown  the  unbiasedness  of  the 
largest-smallest  root  test.  No  study  has  been  made  concerning  and  in 
the  two-sided  case.  Some  progress  in  this  direction  is  made  at  Purdue  University 

A few  other  statistics  arc  available  in  the  literature  which  are  less 
important  than  any  given  above: 

5}  Wilks'  statistic,  Z(s)  - n b.  (1932) 

i«l  1 

0)  Pillai's  harmonic  mean,  ll(s^  « ( £ bT1/*)"1, 

1 i-1  1 

* ( l (l-b.rVs)"1  and  - ( [ f^/s)’1,  (1955) 
i-1  1 ^ i-1  1 

and  7)  Bagai's  statistic,  Y^  « n f.  (1962a). 

i-1  1 
(s) 

Wilks'  statistic  Zv  * , has  been  considered  by  (Gnanadcsikan,  Lauh,  Snyder 
and  Yao,  1965,  Roy,  Gnanadcsikan  and  Srivastava,  1971).  Hsu  (1940a)has  studied 
the  moments  and  distribution  problem  of  Wilks'  statistic.  The  density  of  Zw 
has  been  obtained  by  Sudjana  (1973b)  and  power  comparisons  made  for  the  two-roots 
case  with  those  of  the  four  test  criteria  i)  to  4)  for  each  of  the  three 
hypotheses  I)  to  III),  The  results  show  that  the  test  is  generally  behind  in 
power  to  U<2>.  and  W^  but  has  greater  power  than  the  largest  root.  This 
general  discouraging  aspect  of  the  power  of  the  test  was  confirmed  independently 
by  Hart  and  Money  (1976) . 

Approximations  to  the  densities  of  the  three  harmonic  mean  criteria  have 
been  given  by  Pillai  (1955).  No  adequate  power  studies  have  been  carried  out 
for  the  harmonic  mean  criteria  or  Bagai's  statistic. 
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In  the  tuie-samplc  case  i .>*.  the  c.u.e  ot  a single  covariance  matrix, 

*he  following  hypotheses  are  of  great  interest: 

IV)  . « - (specified)  in  a p-variate  normal  population  N(w,lJ. 

•>  ■» 

V)  In  *“£  ,a*  > 0 unknown  and  £ specified  in  N(u,E). 

For  IV),  Roy  (19S7)  has  suggested  the  largest-smallest  root  test  based  on 
his  union-intersection  principle  of  hypothesis  testing.  Schuurmann  and  Kaikar 
(19**5)  have  studied  the  power  function  of  the  test  in  the  bivariate  case  and 
shown  that  the  test  based  on  the  acceptance  region  (u_i  ij  < ^ < ij)  ■ l-2a 

where  with  s«p)  considered  by  Clemm,  Krishnaiah  and  Kaikar 

(19"2i  enjoys  generally  less  power  than  the  test  studied  by  Thompson  (1963) 
and  llanmara  and  Thompson  (1963)  whose  acceptance  region  is  given  by 
PrU'  « l.  < t < £")■  1-2**,  Prfi,  • l')  ■ 1-a.  However  the  exact  values  of 
fc'  and  t"  were  given  only  for  |*2  and  approximate  t'  and  t"  computed  for  p £ 3 
by  taking  separate  lower  oA  points  for  and  upper  a*  points  for  l so  as  to 
obtain  a 2*  level  test.  Further,  some  exact  powers  have  been  given  by 
Sugiyama  (1972a)  for  the  bivariate  case  based  on  the  largest  root  test  alone 
and  Muirhcad  (1974)  has  obtained  approximate  powers  for  the  same  test  using 
asymptotic  expansions  against  the  alternate  hypothesis  l ■ XI  for  p«2  and  3. 


(Anderson,  1958J  is  a function  of  ire  where  the  g!s  follow  (2.7)  with  s»p  and 
v df.  The  exact  distribution  of  l.  has  been  obtained  by  Nagarscnker  and  Pillai 
(1973)  using  inverse  Mcllin  transform  and  percentage  points  tabulated  for 
p«2(l)10.  The  unbiasedness  of  the  L-tcst  has  been  shown  by  Sugiura  and  Nagao 
(1968).  The  monotonicity  of  power  function  with  respect  to  each  of  the  p ch. 
roots  of  Z Eg1  has  been  established  by  Nagao  (1967)  and  Das  Gupta  (1969). 

For  hypothesis  V),irc  is  a function  of  W * n g./(  f g./p)*5  (Anderson,  1958) 

i«l  x*l 

which  is  the  sphericity  criterion  of  Mauchly  (1940a, b)  who  obtained  the  null 
distribution  for  p*2.  Nagarsenkcr  and  Pillai  (1972,  1973)  have  obtained  the 
exact  null  distribution  of  W in  scries  form  and  tabulations  of  5%  and  1%  points 
for  p*2(l)10  made.  The  unbiasedness  of  the  test  has  been  sh^wn  by  Glescr  (1966), 
Sugiura  and  Nagao  (1968).  Muirhcad  (1976)  has  computed  some  powers  in  the 
two-roots  case  based  on  the  distribution  of  W which  he  has  obtained  as  mixtures 

i 

of  F distributions. 
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In  thv  tUwy  of  jjrmcip.il  component  analysis  (Hotelling,  1933,  193*0 

it  i%  veil  ipifrii  viiat  tae  variances  of  the  principal  components  arc  the  eh. 

r.«u*  *»f  the  covariance  matrix  £.  Since  the  maximum  likelihood  estimates  of 

the  eh.  roots  are  the  eh.  roots  of  the  nte  of  L,  it  is  easy  to  see  that  the 

* 

individual  g^  In  (2.7)  and  its  density  are  of  added  interest.  Thus  the  next 
section  is  devoted  mainly  to  the  study  of  individual  ch.  roots  from  (2.1), 

( 2 . f* j and  (2.7). 


4.  INDIVIDUAL  CHARACTERISTIC  ROOTS 

lour  different  method*  have  been  developed  over  the  year',  for  the  study 
of  the  distributions  of  individual  ch.  roots  as  follows:  a)  Roy-Pillai 

reduction  formulae  b)  Mchta-Krishnainh  pfal.ians  e)  Davis  differential 
equations  (dc'sl  and  Pillai-Pukutomi-Sugiyana  tonal  polynomial  series 

4.1  ROY-PII.LAI  REDUCTION  FORMULAE 


Consider  the  integral  over  the  domain  0 < Xj  <,  x,  £ . . . ^ xk  <_  x £ 1 
of  the  function 


(4.1) 


,r  *xi 


n x?(i-xj  e * a (x.-x.j, 

i*l  11  i!*j  1 3 

where  q,r  •— 1 and  t is  independent  of  the  x's.  Let  q^q,,..^  be  real  nui  rs 
greater  than  -1.  Let  us  denote  by  V(x;qk,...,qt;r;t)  the  psuedo-determir 
(pseudo  in  the  sense  that  the  determinant  has  t#  be  expanded  preserving  the 


order  of  integration) 


(4.2) 


/ xkKfl-xk)re  K dtk.  . ./ c K dxk 


tx. 


tx. 


x«>  \ r tXi  x9  r txi 

/ " x,  (l-x.)r  e 1 dx,.  . ./  x/u-xj  e 1 dx. 
U 1 1 0 


Then  it  is  easy  to  see  that  the  integral  of  (4.1)  can  he  expressed  as  (4.2)  with 
a q+j~l (j*l » • • • >k) . Pillai  (1954,  1956b)  has  proved  the  followi,  g theorem: 
Theorem  4.1.  The  pseudo-determinant 

(4.3)  V{x;qJc,qk_1,...,q1;r;t)  * (q^r+1)"1  (A(k)  + B(k)  + qkC(k)*  tD(k)), 

where  A(k)  « -I0(x;qk;r+l;t)V(x;qkl,...,q1;r;t), 
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As  a special  case  of  Theorem  4.1  (a  second  special  case  will  be  shown 
later  for  the  mgf  of  V^}  putting  t«0,  q^m^j-1, j*l , . . .s,  the  following 
reduction  formula  for  the  cdf  of  the  largest  root  be  is  obtained  (Pillai, 
1954,  1956b) 

(4.4)  (ra+n*s)V(x;m*s-l,...,m;n;0)  « -I0(x;m*s-l;n*l;0)V(x;m*s-2,...,m;n;0) 

s 2 l (-l)!,"',"1I(x;2m*s*j-2;2n+l;0)Vfx;**s-2,iM'j,m+j-2,...,*;n;0) . 

Using  (4.4),  Pillai  (1954,  1956a)  has  obtained  explicit  expressions  for  the 
cdf  of  b^  (s«l,...,s)  in  terms  of  incomplete  beta  functions.  However,  for 
larger  values  of  s computation  is  extremely  prohibitive  and  for  this  reason 
an  approximation  of  the  cdf  of  the  largest  root  at  the  upper  end  has  been 
suggested  by  Pillai  (1954,  1956a,  1965a).  This  approximation  will  be 
discussed  later. 

The  prohlcm  of  obtaining  the  cdf  of  the  largest  root  has  been  investigated 
first  by  Roy  (1945,  1957)  who  gave  explicit  expressions  for  its  cdf  for  s»2,3 
and  4.  Nanda  (194Sa)  also  gave  such  expressions  for  s»2(l)5,  but  no  general 
expression.  Pillai  and  Dotson  (1969)  has  proved  the  following  lemma: 

th  th 

lemma  4.1.  If  bj  and  b are  the  i 1 and  (s-i+1)  roots  (i»l,...,s) 
where  bj,...,b  follow  the  density  (2.1),  then 

(4.5)  Pr(l>.  < x;  m,n)  « 1-Pr(bs_i+1  <,  l-x;n,m), 

where  on  the  right  hand  side  of  (4.5)  m and  n arc  interchanged. 

Roy  (1945)  has  also  given  an  expression  for  the  cdf  of  the  i^  root 
(i«l,...,s).  However  the  expression  is  not  correct  and  even  the  total 
probability  for  the  whole  range  docs  not  equal  unity.  Pillai  and  Dotson 
(1969)  have  given  the  co-r.ct  reduction  formulae  for  the  smallest  root, 
second  largest  and  soeo.  i smallest  roots.  (Formula  (4.5)  is  not  useful  when 
computations  for  larger  valuus  of  m cannot  be  easily  made  as  is  usually  the 
case.)  They  have  also  given  upper  5*  and  1%,  points  for  the  smallest  root  for 
p *2  and  3 and  middle  root  for  p<*3. 


Returning  to  the  pro!  l>vi  of  the  extreme  roots,  Pillai  (1965a)  has  obtained 
general  expre  . uppr.Kt*  it  ing  ».*  the  upper  end  the  cdf  ot’  the  largest  of 
s turn -null  eh.  rootN.  these  expression**  .ire  the  on  tv  source  for  computing  the 
ulf  of  the  largest  eh.  root  for  larger  values  of  s,  and  are  given  below, 
ill  s even,  (he  approximation  to  1’rfh,  ^ xi  for  even  values  of  s is  given  by 

i 

I I.O)  1 + I ( -1 ) l Ifl|x;ms.i;n+1;0),  kmfs«0,  wherc 

i*  1 


(m*n«-s-i+l)"ItC(s,mln}V(l;m+-s-l , . . . ,m+s-i+l ,m*s-i-l ,. . . , 


m;n;0)-(mts-Ul)kmfs_.  + 1), 


C(s-l,ra,n  )V(l;m*s-l,...,n*s- 


(2w»s-j»i: 


x*l,m+s-x-l,...  ,m;n;0)  « 


(Li)  s odd.  The  approximation  to  Pr(bs  £ x)  for  odd  values  of  s is  given  by 

<4-7>  H&rtSifi * 

MHt  1 J 

The  approximation  neglects  terms  involving  (1-x)  , r>l,  in  the  development 

of  the  cdf  obtained  by  repeated  application  of  (4.3).  Pillai  (1954,  1956a,  1957, 
19t>0)  gave  the  upper  5 and  1%  points  of  bs  for  s*2  and  5 using  his 
approximation  formulae.  Sen  (1957)  computed  similar  upper  percentage  points 
for  three  roots  and  Ventura  (1957)  for  four,  both  following  Pillai's  method. 
Pillai  and  Bantcgui  (1959)  gave  such  tabulation  for  s»6.  All  these  percentage 
points  were  given  for  values  of  m*0(l)4  and  n varying  from  5 to  1000.  Further 
Jacildo  (1959)  extended  the  tables  for  s»2  and  s»3  for  values  of  m*5,7,10  and 
15  and  the  same  range  of  values  of  n as  before, Pillai  (1960)  has  published  all 
these  percentage  points  for  s*2(l)6.  Further, Pillai  (1964a)  has  given  such 
percentage  points  for  s»7  for  values  of  m«0(l)5,7,10  and  values  of  n as  before. 
For  all  the  computations  up  to  s*7  the  approximation  was  obtained  for  each  s 
since  the  general  expressions  (4.7)  and  (4.8)  were  obtained  only  in  1965  when 
percentage  points  were  also  given  for  s*8,9  and  10  for  values  of  m*0(l)5,7,10,15 
and  n as  before  (1965a).  Similar  percentage  points  for  s*ll  and  12  arc  given 
in  Pillai,  (1970)  and  for  s*14(2)20  in  (1967b).  Further,  percentage  points  for 
s*2(l)10  and  s=13(l)20  are  available  in  two  mimeograph  reports  in  the  Department 
of  Statistics  at  Purdue  University  (Pillai,  1966b,  No. 76,  1966a,  No. 72). 

Nanda  (1951)  has  given  upper  5 and  1%  points  of  the  largest  root  for  s*2 
and  very  small  values  of  m and  n(m=0($)2;n=£(£)10.  FosteT  and  Rees  (1957)  have 
tabulated  the  percentage  points  (80,85,90,95  and  99)  of  the  largest  root  for 


s*2  m •*-<*. 5,  0(l)il  ami  n«l(  1)13(5)49,59,79.  Foster  (1957,  1958)  has  further 
extended  these  tables  for  values  of  s*3  ami  4.  The  arguments  they  have  used 
for  tabulation  arc  the  degrees  of  freedom.  Heck  fl960)  has  given  some  charts 
of  upper  5,  2.5  and  1}  points  for  s*2(l)5,  m»-f,  0(1)10  and  n > 5. 

In  the  one  sample  case  reduction  formulae  for  individual  roots  could  be 
obtained  from  these  discussed  above  by  the  method  used  for  obtaining  (2.7) 
from  (2.1).  For  example,  for  the  largest  root  the  reduction  formula  for  the 
one  sample  case  can  be  obtained  from  (4.5)  in  the  following  manner: 

Let  W(y;qk, . < . ; -1)  denote  the  pseudo-determinant  obtained  from  (4.2) 
by  putting  r*0,  t*-l  and  x^  changed  to  y^(i*l, . . . ,k)  and  x to  y where 

0 * >’i  *.•••*  yk < *•#  t^en 

(4.8)  h (y ;qk , . , . ,q^ ; -1)  * “iQCyjq^I  *l)h()'jqj._j , • • • ,qj  J“l) 

* 2 J (~i)  ^ i »“2) K(y»qj(_ j , • • « ,qj+ j »qj_ j » • • • ♦tij » "1) 

♦ qk  W(y;qjjl,qk_l,...,ql;-l), 

y ty  ty  y 

where  I(y;q,t)  * / yj  e 1 dyt  and  IQ(y;q;t)  - yj  e *|  . 


Hence  the  cdf  of  g can  be  obtained  from  that  of  b„  and  similarly  those  of  other 
individual  roots.  Nanda  (1948b)  has  given  explicit  expressions  for  the  cdf  of 
the  largest  root  g for  s*2(l)5  and  similarly  the  smallest  root  and  middle  root 
for  s*3.  However  since  the  computation  is  prohibitive  for  larger  values  of  s 
Pillai  and  Chang  (1968,  1970)  and  independently  Hanumara  and  Thompson  (1968) 
obtained  an  approximation  to  the  cdf  of  gg  (approximating  at  the  upper  end) 
starting  from  Pill ai * s approximation  for  the  cdf  of  bg  in  (4.6)  and  (4.7)  and 
using  the  method  for  deriving  (2.7)  from  (2.1).  The  approximation  is  as  follows: 

(i)  s even.  The  approximation  of  Pr(gs  <.  y)  for  even  s is  given  by 

f4.9)  ! . y"*'-1  e-y,  Cs  - 0 

1*1 

whore  Cs-i  * .Hj  (»*-m-(w-w)C,.w 

(ii)  s odd.  The  approximation  to  Pr(gs  <. y)  for  odd  s is  given  by 

«•>»>  Tc-tf  c.i  ^ °-y- 

Pillai  and  Chang  have  tabulated  (1968,  1970)  upper  10,5,2.5,1  and  0.5%  points 
of  gs  for  s*2(l)20  and  df  ranging  from  2 to  200.  Hanumara  and  Thompson  (1968) 


have  tubulated  upper  percentage  points  as  above  except  1G*  for  s»2(l)10  and 
also  lower  percentage  points  for  the  smallest  root  derived  in  view  of  lcwna 
('1.1} . burlier,  Thompson  (1962}  had  tahulatcd  upper  percentage  points  for 
s«2. 

4.2  MEIITA-KRISIINAIAH  P FAFF I AN  METHOD 


Mehta  (1960,  1967}  has  developed  an  alternate  method  of  reduction  of  the 
pseuedo-determinant  in  terms  of  double  integral*  which  was  extensively  used 
with  extensions  by  Krishnaiah  and  associates  (Krishnuiah  and  Chang,  1970,  1971b, 
Krishnaiah  and  Kaikar,  1970,  1971b}.  The  reduction  is  as  follows:  Let 


(4.11) 


A 

»(x 1 ,x";q. , . . . ;^)  * /..../  *•  • • H C'S'Cx. )dx*. } , 


M * 

where  ^o^l  * * * * [’^1 1 * * * * I (xj  )i»  j*l » • • jk|  • Further^lct 

x"  u 

Fj(x',x")  * / F (x',u)u  #(u)du,  F (x* ,u)  « / vTy(v)dv. 
r x'  r r x' 

and  f*(x',x")  - fJ(x,,x")-fJ(x’,x").  In  addition,  let 

ncx',x»';q1,..,q2f;«)  - | (f^(x*  ,x"))i,j-l 2f|* 

and  (x1  ,x  ' ;q^ , . . . ,q2£+j JV1)  * 1 (fq^  (x*  ,x"))i  p j*l , . • . ,t-l  ,t*j  t , j2f+l  | ■ . 

Note  that  both  l)  and  G arc  pfaffians.  (If  T is  a skew  symmetric  matrix  of 
J * 

even  ordcr^  |T| is  called  a pfaffian.) 

Lemma  4.2.  Let  ^(x)  be  a function  such  that  the  integral  given  in  (4.11) 
exists  and  let  x'  < x"  be  real  constants.  Then 

(4.12)  V(x,,x";q1,...,qIc;V)  * D(x,,x";q1,...,q2£;»J0/f!  when  k«2f 

2f+l  . 

(4.13)  * l (-If  JFq  (x',x")G.(x',x»;q1,...,q2f+1;<»)  when  k-2f*l. 

Mehta  (1960)  proved  (4.12)  in  the  special  case  when  q.*i-l,  i»l,...,k 
_xn  r 1 

and  ^ (x)  =«  e x . Krishnaiah  and  Chang  (1970  , 1971b)  proved  the  above 

lemma  using  Mehta's  method  (1960)  for  the  case  when  q^*i-l ,i*l , . . . ,q.  Lemma 

(4.2)  is  due  to  Krishnaiah  and  Waikar  (1970,  1971b).  It  is  obvious  that 

(4.11)  multiplied  by  appropriate  coefficients  include  as  special  cases  the 

cdf's  of  the  largest  root  from  (2.1)  and(2.7). 

Now  in  regard  to  the  root,  x^,  i ^ j < k, 

(4.14)  Pr(Xj  £ x)  = Pr(Xj+1  < x)  + Pr(Xj  < ...  < x^  < x < x^  < ...  < x^.) 


which  gives  a recurrence  formula  for  the  cdf  of  intermediate  roots.  Once  the 
cdf  of  the  largest  root  is  known  it  is  enough  to  evaluate  the  second  term  on 
the  right  side  of  (4.14).  (A  similar  approach  knowing  the  cdf  of  the  smallest 
root  first  is  also  available.)  Now  developing  Vq(  ) in  (4.11)  by  the  first  j 
columns  (noting  that  in  the  densities  (2.1)  and  (2.7)  q^«i-l)  and  using  lewaui 
4.2  to  evaluate  the  integral 

(4.15)  Pr(Xj  * x)  ■ Pr(Xj4l  < x)  ♦ c Ej  ± Wa,x;qj , . . . ,qJ.;W  aCx.bjqJ,^ , . . . ,q£;^» . 

where  (a,h)  is  the  rang.?  of  the  variables  and  c the  constant  in  the  joint  do  my 
in  question  and  qj  * ...  * qj.  is  a subset  of  the  integers  0,1,..., k-1  and 

^r+1  v •••  ' ^k  C^e  C0!nP*t,,r‘cnt  set>  *•]  denotes  summation  over  all(^)  possible 

choices  of  qj  * ...  < qj,.  Further,  Krishnaiah  and  Kaikar  (1970,  1971b)  have 
obtained  the  joint  density  of  any  suhset  of  consecutive  ordered  roots  in  a 
similar  manner. 

Based  on  the  joint  distribution  of  the  largest  and  smallest  roots  bj  and  bs 
obtained  by  Krishnaiah  and  Chang  (1970,  1971b)  by  methods  described  above, 
Schuurmann  Kaikar  and  Krishnaiah  (1973)  have  obtained  values  of  x such  that 
Pr(l-x  < b,  < b < x)  * l-o  for  o * 0.10,  0.05,  .025,  0.01;  »»0 (1)5,7,10,15; 
n»5(l) 10(2)20(5)50  and  s»2(l)10.  Since  optimal  choice  of  x'  and  x"  is  not 
known,  the  authors  carried  out  these  tabulations  fer  facilitating  Roy's 
two-sided  test  of  hypothesis  I) , Tables  for  the  upper  percentage  points  of 
b^  has  also  been  given  by  Schuurmann  andKaikar  (1972) . Further  Clemm, 

Krishnaiah  and  Kaikar  (1972)  have  tabulated  the  values  of  u such  that 
Pr(u_1  £,is<.u)  ■ 1-a  where  * 2g^(i»l, . ..  ,s)  for  values  of  a ■ 

0.05,  0.025,  0.01,  0.005,  s«2(l) 10(2)20,  and  df  « (p*l) (1)20(2)30(5)50. 

This  again  has  been  to  facilitate  Roy's  two-sid**d  test  of  hypothesis  IV)* 

Thompson  (1962)  had  computed  the  same  for  s»2.  They  have  also  provided 
various  lower  percentage  values  u,f  Further,  upper  percentage  points  of 
the  individual  roots  ■ 2g.(i«l,...,s-l  for  s«2(l)10  and  the  above  values  of 
df  have  been  given  by  Clemm,  Chattopadhyay  and  Krishnaiah  (1972). 


4.3  DAFIS  DE  METHOD 


A third  approach  is  due  to  Davis  (1972a,  1972b)  who  has  shown  that  the 
marginal  density  functions  of  bjs(gjs)  satisfy  a system  of  ordinary  differential 
equations  of  Fuchsian  type.  Let  ^,v2  >,  p in  (2.1)  so  that  s*p.  For  notational 
convenience  let  c.*l-hj(i*l,...,p) . Further  define  DT(t,c)  * {0  < xr  < ...  < xfc 
<c<xtj<...<x^<l).  Then  the  marginal  density  function  f^Cc)  of  ct  is 
given  by 


I i $niV,,v.(X.,...,x  .»d,X,,..,x  »]dx 

1 V'1(t.c)  r *•  1 1 C_1  c p 1 - 


P-1 

where  dx  ■ n U , It  is  proportioned  to 
* i*l  xi 

f(v,-p-l)  ^lv.-p-U  P-1 

c “ tl-cj  J *(x)  n (e-x.Jdx, 

Pp_1lt,e)  * i»l  1 * 

where  *>  denotes  ♦ , , , , a»d  $ denotes  the  joint  density  of  the 

p-tiiyijVj-i  P,V2|V1 


c!s.  Again,  define 

yc,x)  « «(x)Sa(c*xci(u)...Cc-xo(p_1-r))/r-0,l p-1), 

the  summation  being  extended  over  the  (p“*)  selections  of  integers 

a(l)  * ...  « a(p-l-r)  from  che  set  l,2,..,p-l.  When  r ■ p-1,  the  sum  is  taken 

to  be  unity.  Now  introduce  the  p functions 

K r(c)  * / n i M«»*)dx,  (r»0,l,...p-l). 
t,r  0P  (t,c)  r 

The  system  of  differential  equations  derived  by  Davis  is  the  following: 

M.lbJ  c(l-e) l^r  - f.(p-r)(v1*v2-p*r-l)Lt,r_1 

♦ i r[(l-c){v,-p+r)-e(v.-p*r)].l, 

••  * W|l 

♦ $(r+l)(r*2)c(l-c)lt  ^ t (r«0,l p-1), 

where  L £i.  X 0.  Further,  it  is  convenient  to  introduce  H «(l-c)rL 

1 1 • i p r 1 1 r 

(r«0, . . . ,p-l)  and  to  express  (4.16)  as  a matrix  de  for  ii  ■ (IIq,...,H  j) 
given  by 

(4.17)  dH/dc  - (c”1A*(l-c)”^C]H, 


where  matrices  A and  C arc  given  in  (1972a).  The  de  (4.17)  is  of  the  Fuchsian 

•m  •• 

type  with  regular  singularities  at  c=0,l  and  ».  The  marginal  density  of  the 
largest  root  has  been  solved  as  a power  scries  which  coincides  with  that 
obtained  by  Pillai  (1967a)  and  independently  by  Fukutomi  and  Sugiyama  (1967). 
Pillai's.  approximation  ;o  the  distributions  of  the  largest  (smallest)  root  in 
(4.6)  and  (4.7)  arc  interpreted  as  exact  solutions  the  contributions  of  higher 
order  solutions  being  neglected.  No  approximations  of  a useful  nature  have 
been  suggested  for  the  distributions  of  the  other  ch.  roots.  Further,  Davis 
(1972b)  has  given  explicit  expressions  of  the  marginal  distributions  of 
individual  roots  bjs  for  p £ 5.  A de  is  applied  recursively  to  construct  the 
distributions  for  p from  those  obtained  from  p-1.  Results  also  have  been  derived 
for  = (2g!s) . The  results  for  fjs  in  (2.6)  would  follow  from  those  of  b!s. 
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Finally  the  system  of  lie's  discussed  above  is  related  to  Davis'  de's  for 
Hotelling's  trace  (Davis,  1968,  1970a,  1970c)  ar.d  those  for  Pillai's  trace 
(1970b). 

The  work  of  Davis  and  earlier  authors  have  been  generalized  in  various 
directions  by  Eckert  (197S)  where  the  >y  function  in  (4.11)  includes  the  beta, 
gamma,  normal  and  in  addition,  a class  of  V'-functions  given  by  weight  functions 
for  classical  orthogonal  polynomials  to  include  Pearson  system  of  curves 
distinguishing  three  more  types.  Ilis  reduction  formulae  generalize  the  work 
of  Roy  (1945,  1957) ,Nanda  (1948a, bJ^Pillai  and  Dotson  (1969)  and  Pillai 
(1954,  1956b)  when  t»0.  He  also  generalizes  the  work  of  Davis  (1972a,  1972b) 
in  the  context  of  the  Pearson  system.  Further  Eckert  defines  an  operator  to 
find  ail  the  cdf's  of  the  individual  roots,  and  in  regard  to  the  normal 
finds  the  basic  algebraic  form  for  the  distribution  of  the  largest  root  in 
terms  of  a certain  space  of  functions  and  one  can  pass  on  from  the  distribution 
of  the  largest  root  to  those  of  the  other  roots  through  certain  operators. 

4.4  P-F-S  ZONAL  POLYNOMIAL  SERIES 

A fourth  approach  which  has  not  been  very  useful  in  the  null -case  in  view 
of  the  better  approaches  above  but  quite  useful  in  the  non-null  case  cf  the 
distribution  problem,  is  the  scries  approach  involving  zonal  polynomials 
(James,  1961,  1964,  llua,  1959,  Subrahmaniam,  1974,  see  Section  9)  which  are 
homogeneous  symmetric  polynomials  in  the  ch.  roots  of  a real  symmetric 
matrix  (Hcrmitian  matrix  in  the  complex  case).  The  density  of  the  largest 
root  bs  has  been  obtained  by  Pillai  (1967b)  and  independently  by  Fukutomi  and 
Sugiyama  (1967)  as  infinite  scries  of  zonal  polynomials.  The  density  of  the 
smallest  root  follows  from  Lemma  (4.1).  Sugiyama  (1967b)  further  gives  the  cdf 
of  bs(bj)as  a series  in  powers  of  x((l-x)).  The  distributions  of  f^ffj)  follow 
by  transformation  used  to  obtain  (2.6)  from  (2.1).  For  the  one  sample  case, 
Sugiyama  (1966,  1967a)  has  obtained  the  distribution  of  Jts«2gs  ant*  t*'e 
density  of  the  vector  corresponding  to  useful  in  principal  component 
analysis  using  some  results  of  Tumura  (1965),  and  has  given  an  approximation 
(1972b)  to  the  distribution  of  4 from  the  zonal  polynomial  view  point. 

Further,  Krishnaiah  and  Chang  (1971a)  have  obtained  the  distribution  of  the 
smallest  root  in  terms  of  a finite  series  of  zonal  polynomials  when  m is 
a non-negative  integer. 
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Till.  TWO  TKA13.S 


In  this  sivtion,  tlu*  two  tract's,  those  ol  Hotel  1 inn's  and  I'illai's, 
will  he  considered  in  some  detail.  l;or  the  study  of  the  moment,  PI  Hal 
(1954,  1956b)  has  obtained  a recurrence  relation  concerning  the  mgf  of 
which  is  a special  case  of  (4.3)  in  Theorem  4.1  with  x*l  and  q^  ■ 
ra+j-1 ()«1, . . ,sj  as  given  below. 

*v(S)  fsl 

(5.1)  (m+n+s-t)F.(ccv  )+tl:(VlSjecv  )-[2C(s,m,n)/C(s-2,ra,n|) 


Sf  l-l)s“j'1lCl;2rafr.+j-2;Jn+l;2t)E[(Eb.,...,b  ),  ,ctV 
j«l  1 s“3  1 


(s-2) 


). 


t*  fh 

where ^bj .. .bs  ^ 1 denotes  the  (s-j-1)  esf  in  (s-2)  variables 

Pillai  studied  the  first  four  moments  of  using  (5.1)  and  suggested  a beta 

function  approximation  to  the  null  distribution  of  (Pillai,  1954,  1955, 

1957,  1960).  Pillai  and  Mijares  (1959)  gave  an  alternate  method  of  expressing 

the  moments  of  and  then  the  general  expression  for  the  fourth  moment 

(Ting,  1959,  Pillai,  1960)  which  was  obtained  by  Pillai  only  for  s«2,3  and  4 

earlier.  The  first  four  moments  were  used  to  obtain  approximate  upper  5% 

and  1*»  points  of  (Pillai,  1957,  1960)  for  s»2(l)8  and  various  values  of 

m and  n.  (The  method  of  using  moment  quotients  to  obtain  percentage  points 

is  available  in  Pearson  and  Hartley  (1956)).  Further, Mijares  (1964a) has 

tabulated  the  upper  and  lower  5%  and  1%  approximate  percentage  points  of 
Ts ) 

Vv  J for  s«2(l)S0  using  the  moment  quotients. 

In  regard  to  Hotelling's  trace,  as  mentioned  in  the  previous  section, 

Pillai  (1954,  1956b)  has  obtained  a theorem  similar  to  Theorem  4.1  concerning 

the  pseudo-determinant  for  (2.6)  in  terms  of  the  £!s  from  which  a recurrence 

(s)  1 

relation  for  the  Laplace  transform  of  lr  * corresponding  to  (5.1)  follows  as  a 

(s') 

special  case.  Pillai  (1954,  1956b)  has  obtained  the  first  four  moments  of  U 

and  suggested  an  F approximation  to  the  null  distribution  of  the  statistic. 

Pillai  and  Samson  (1959)  have  obtained  approximate  upper  5%  and  10%  points 

for  s=2,3  and  4 using  the  moment  quotients  and  Pillai  (1957,  1960)  has  given 

such  approximate  percentage  points  for  s*2 (1)8  and  various  values  of  m and  n. 

Further,  Pillai  (1964b)  has  obtained  the  following  lemma  which  enables  one  to 

obtain  the  moments  of  from  that  of  V^(and  vice  versa). 

Lemma  5.1.  Let  and  U.^  denote  the  ith  elementary  symmetric  functions 

i,m,n  i,m,n  .. 

in  the  s g|s  in  (2.1)  and  s f!s  in  (2.6)  respectively.  Then  the  kin  moment 

u'{U^  } is  derivable  from  u,!.{V?s^  } by  making  the  following  changes  in  the 

k ljiUjH  k i j m j n 
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expressions  for  the  Utter:  (a)  Multiply  by  -1  all  terms  except  the  term  in  n 

in  each  linear  factor  involving  n and  (b)  change  n to  m*n*s*l  after  performing  (a). 

Other  approximations  to  the  distribution  of  Hotelling's  trace  generally 
based  on  th'/  moments  have  been  given  by  Pillai  and  Young  (1971)  Tlku  (1971)  and 
Hughes  and  Saw  (1972). 

In  regard  to  the  exact  distributions  of  the  traces,  there  has  been 

considerable  work  but  no  useful  general  forms  have  yet  been  obtained.  Nanda 

(19S0)  has  derived  the  cdf  of  in  the  special  case  of  m>0  and  s«2  and  3. 

The  method  was  to  study  the  mgf  and  carry  out  inversion  involving  the 

distributions  of  the  weighted  sums  of  two  (or  more)  independent  random 

variables  of  the  form  (2n»ki*l)(l-yi)“n+lci,0  < <1.  Pillai  and 

Jayachandran  (1970)  have  extended  the  method  farther  using  (4.8)  for  reduction 

of  pseudo-determinants  arising  in  the  derivation  of  the  mgf  and  obtaining 

explicit  expressions  of  the  cdf  of  for  s-3  and  integral  values  of  m < 3 

fsl  “ 

and  s>4  and  m«0  and  1.  Further,  exact  upper  percentage  points  of  Vv  have 

been  computed  for  values  of  s and  m given  above  and  selected  values  of  n. 

PI 

Earlier,  Mikhail  (1965)  has  obtained  the  exact  density  function  of  . 

In  regard  to  Hotelling's  trace,  U^,  the  distribution  for  s>2  has  been 
obtained  by  Hsu  (1940a)  which  has  the  same  form  as  given  by  Hotelling  (1951), 
in  terns  of  an  incomplete  beta  function.  The  density  can  be  expressed  in  the 
form  of  a hypergeometric  series.  Pillai  and  Young  (1971)  and  Pillai  and 
Sudjana  (1974)  have  developed  the  Laplace  transform  of  U(s)  as  a 
psuedo-deterninant,  and  using  appropriate  reduction  formulae  and  inversions 
involving  convolutions  of  two  (or  more)  independent  random  variables  (special 
type  of  F variables),  obtained  the  exact  distribution  of  for  p*3  and  4 
and  small  nen-negative  values  of  m. 

Davis  (1968)  has  shown  that  the  density  function  of  in  the  null  case 
satisfies  an  ordinary  linear  differential  equation  of  Fuchsian  type  and  has 
computed  percentage  points  (Davis,  1970c)  by  analytic  continuation  of 
Constantine's  series  distribution  of  (Constantine,  1966)  which  is 
convergent  only  for  |U^|  < 1.  Actually  Constantine's  series  in  the  null 
case  reduces  to  the  relevant  solution  of  the  Je  ir  the  unit  circle  about 
U^«0.  The  tabulation  is  made  of  5\  and  1%  points  for  v2U^sVvj  for  s=3  and 
4 and  various  values  of  and  v2<  Pillai 's  approximate  percentage  points 
(1957,  1960)  have  been  generally  shown  to  have  three  decimal  accuracy  except 
for  small  values  of  v„. 


in  the  null  case 
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Further,  Davis  (1970b)  has  shown  that  satisfies  a hide  of  Fuchsian 
type  which  is  related  by  the  following  simple  transformation  to  his  de  for 
U^).  u^**-V^,  Vj-'S-Vj-Vj+I.  The  transformation  brings  out  the  relationship 

between  the  moments  of  and  (See  Lemma  5.1  of  Pillai).  Again,  Davis 
has  found  that  Pillai's  approximate  percentage  points  of  (Pillai,  1957, 
I960)  are  accurate  to  four  decimal  places  except  when  and  v2  are  both  small. 
Davis  has  also  shown  that  the  density  f (V^)  ■ f (s-V^)  which  can 

Vv2  v2»vl 


be  used  to  obtain  lower  percentage  points. 

Again,  Krishnaiah  and  Chang  (1972)  have  studied  the  exact  distributions 


of  and  using  the  pfaffian  method  discussed  in  the  previous  section 

to  obtain  the  Laplace  transforms  as  linear  combinations  of  the  products  of 

certain  double  integrals.  Further,  for  the  double  integrals  have  been 

developed  by  Schuurmarat,  Krishnaiah  and  Chattopadhyay  (1973)  as  linear 

combinations  of  incomplete  gamma  functions  and  inverted  to  obtain  the  density 

of  explicitly  for  s»2(l)6  and  several  values  of  v.  and  v,  The  authors 

fsl  1 

have  given  the  exact  upper  percentage  points  for  Vv  1 for  s»2(l)S,  o*. 01,. 025, 

.05  and  .10  for  m and  n»0(l) 10(5)25.  They  have  confirmed  the  comments  of  Davis 

concerning  Pillai's  approximate  percentage  points.  As  regards  Krishnaiah 

and  Chang  (1972)  have  considered  explicitly  only  the  case  of  s»2  and  v^S. 

In  regard  to  asymptotic  distributions  of  and  asymptotic 

2 fs  ^ 

expansions  for  Hotelling's  TQ(»v2Uv  ')  both  for  the  cdf  and  percentile  have 

been  given  by  Ito  (1956)  as  a chi-squared  series  up  to  order  the  first  term 
•)  ^ 
being  . An  independent  derivation  of  Ito's  expansion  is  given  by  Davis 

(1968,  19)oa).  Davis  (1971)has  extended  Ito's  expansions  to  order  vl3  . 


,(*) 


Muirhead  (1970b) has  considered  a general  asymptotic  expansion  for  functions 

satisfying  a certain  system  of  partial  de's,  of  which  the  distributions  of 

y(s),  y(s)  an(j  iargCSt  root  in  the  one  sample  case  fall  as  special  cases. 

_2 

Muirhead  explicitly  obtains  asymptotic  expansions  up  to  order  v2  for  the 
cdf's  of  and  the  largest  root.  Further,  Davis  (1970b)  has  obtained  an 
expansion  for  up  to  order  v23. 

Pillai  (Pillai,  1973,  Pillai  and  Sudjana,  1974)  has  suggested  an  F-typc 
series  form  for  which  is  exact  for  s*2  in  the  null  case.  The  exactness 

has  not  been  verified  beyond  s*2.  The  scries  involves  zonal  polynomials. 
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6.  WILKS'  CKUhKIOM 

Let  A.  and  A-,  be  positive  definite  symmetric  matrices  of  order  p,  having 
»»  * ■»  “• 

independent  Wishart  distributions  W(p,vj ,E) , i*l,2,  respectively,  (Anderson, 
1958,  Kshlrsagar,  1972,  Roy,  1957).  Let  A%CBC'  and  A.*A,»CC',  where  C is  a 
lower  triangular  matrix.  Then  it  has  been  shown  (Hsu,  1939)  that  B and  C 

•m  mm 

are  indcpdcntly  distributed  and  the  density  function  of  B is  given  by 

(Khatri,  1959,  Olkin  and  Rubin,  1964). 

. $(v2-p-l) 

(6.1)  Cl(p,vi,v2)|B|  1 | 1-Bj  2 

where  c^p.v^v,)  . n"P(P‘lJ/4  n^r^^-itii/r^Vj.itDjr^.iti)]). 

The  density  (6.1)  is  known  as  the  multivariate  beta  distribution.  If  B is 
further  diagonalized  by  an  orthogonal  transformation  anr!  integration  is 
carried  out  with  respect  to  the  random  variables  in  the  orthogonal  matrix, 
we  obtain  the  joint  density  of  the  roots  b!s  in  (2.1).  The  method  can  be 
modified  suitably  for  the  case  of  A^  being  at  least  positive  semi*definite. 
(Sec  Section  9,  cq  (9.10)). 

Again  let  L - I-B.  Then  Wilks'  criterion  W(p)  ■ |l|.  Wilks  (1932)has 
•<•  «•  •» 

shown  through  the  uniqueness  of  the  moments  easily  derivable  from  (2.1)  that 
the  density  of  W^  could  be  written  as  a product  of  p independent  beta 
variables.  (Sec  Theorems  8.5.1  and  8.5.2  of  Anderson,  1958).  Now  let  L«TT', 


th 


where  T is  a lower  triangular  matrix.  Then  |l|  » n t..,  where  t..  is  the  i 

i-1  11  11 

diagonal  element  of  T.  Kshirsagar  (1961)  showed  that  t^(i»l,..,p)  are 

. 2 

independently  distributed  and  that  t^  follows  the  distribution 

(6.2)  f.(t:.)  « (t2.)  2 (1-tT.)  1 /B^(v2-i+l),^v1),0<t2.  < 1. 

Further,  Pillai  (1964c,  1966c) has  shown  that  the  p(p+l)/2  distinct 
elements  of  L can  be  transformed  to  p(p*l)/2  independent  beta  variables  and 
since  * tr  B and  * tr(I-B) "^-p,  an  alternate  method  of  computing 
the  moments  of  these  statistics  has  been  given  using  the  independent  beta 
variables.  Again,  Khatri  and  Pillai  (1965)  has  given  a method  of  obtaining 
the  densities  of  the  independent  beta  variables  in  successive  stages. 

In  regard  to  Wilks'  criterion,  Wilks  (1935)  has  obtained  the  exact  null 
distribution  in  the  form  of  a (p-l)-fold  multiple  integral  which  he  evaluated 
explicitly  for  p»l,2,  for  p=3  and  v^=3,4,  and  p»4  and  v^*4.  Further, 
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asymptotic  approximations  have  been  considered  by  several  authors.  For  example, 
Bartlett  (1938,  1947a)  obtained  a chi-squared  approximation  to  -C  log  W^whcre 
C ■ v,  - f(p-v,+l).  Wald  and  Brookner  (1941)  developed  an  asymptotic  expansion 
which  was  further  modified  by  Rao  (1948)  to  obtain  the  first  three  terms  of  a 
more  rapidly  convergent  series  for  the  cdf  of  -C  log  W^.  Further,  Box  (1949) 
gave  asymptotic  approximations  to  functions  of  general  Ir  statistics  which 
include  Kao's  as  a special  case.  Rut,  in  addition,  Rao  (1951,  1952)  has  given 
a second  approximation  as  a beta  series. 

Schatzoff  ( i 966b)  has  given  a method  for  obtaining  the  exact  distribution 

of  using  the  representation  of  -log  as  a sum  of  independent  variables 

and  taking  successive  convolutions.  However,  he  has  not  given  explicit 

expressions  for  the  density  or  cdf  of  W^  while  Pillai  and  Gupta  (1969)  have 

given  both  for  values  of  p up  to  6.  Schatzoff  used  Theorem  8.5.1  of 

Anderson  while  Pillai  and  Gupta  used  Theorem  8.5.2  (Anderson,  1958). 

■> 

Schatzoff  has  tabulated  the  factors  for  converting  x“  percentiles  to  exact 

P'S 


percentiles  of  -C  log  W^  for  p»3(l)8  and  values  of  v.  such  that  pv,  <,  70. 
Pillai  and  Gupta  (1969)  have  lifted  this  restriction.  Unlike  Consul  (1966c) 
who,  following  Nair  (1938),  has  used  inverse  Mcllin  transform  to  obtain  the 
distribution  of  for  p up  to  4 as  infinite  series,  Pillai  and  Gupta  (1969) 
have  given  the  same  in  finite  series  form  except  when  both  p and  are  odd  in 
which  case  the  series  is  infinite.  Further  for  p»3(l)6,  they  have  extended 
Schatzoff* s tables  for  selected  values  of  Vj  >,  11  and  <,  22. 

Lee  (1972)  has  studied  the  density  and  cdf  of  W^  when  p or  Vj  is  even 
with  a view  to  simplifying  further  the  expressions  obtained  by  earlier 
authors  for  numerical  computation.  When  p and  are  both  odd,  the  expressions 
given  are  in  terms  of  simple  integrals.  Asymptotic  expansions  also  have  been 
developed.  Tables  of  the  chi-squared  correction  factors  have  been  extended 
to  cover,  in  conjunction  with  previous  tables,  values  of  p £ <.  20  and 

pvj  < 144  with  omission  when  p or  is  odd  and  > 10. 


M;i  thii  i (lU71a,  b)  has  derived  the  null  distribution  using  inverse  Moil  in 
transform,  theory  of  residue  and  psi  and  seta  functions  and  computed  upper  5* 
and  1«  points  for  for  selected  values  of  the  arguments.  Muthui  and  Hathie 
(1971b)  have  obtained  the  exact  null  distribution  of  in  the  general  case 
in  terms  of  simple  algebraic  functions  which  can  be  computed  without  much 
difficulty.  The  density  and  the  cdf  arc  both  given  and  some  special  cases  are 
obtained  from  the  general  form  which  coincides  with  the  results  of  Consul. 

The  method  of  derivation  again  is  taking  the  inverse  Mcllin  transform.  But 
before  taking  the  transform  of  the  ratios  of  the  gamma  products  in  the  kth 
moment  of  elimination  of  the  gammas  by  cancellation  of  common  factors 

has  been  made  and  splitting  the  factors  with  the  help  of  partial  fraction 
methods  has  been  carried  out.  Results  are  given  separately  for  the  four 
cases  of  p and  taking  odd  and  even  values.  Further,  Bagai  (1972a)  has 
alt*o  expressed  the  null  distribution  of  W^  in  terms  of  Meijer's  C-function 
using  inverse  Mcllin  transform. 

7.  OTHER  STATISTICS 

In  this  section,  the  three  statistics  5)  to  7}  of  Section  3,  name'y, 

Wilks'  statistic,  Pillai's  harmonic  mean  criteria  and  Bagai's  statistic  will 
be  considered  first  and  then  statistics  for  the  onc-samplc  case. 

Wilks'  statistic,  Z^s\  considered  by  Hsu  (1940a)  and  Law Icy  (1956)  has 
been  studied  by  Bagai  (1964b)  and  Consul  (1964b)  obtaining  the  limiting 
distribution  using  inverse  Mcllin  transform  for  s»2(l)8.  Ito  (1962)  has 
discussed  Z(s)  as  also  Troskic  (1966,  1972)  and  dc  Waal  (1968,  1970).  Pillai 
and  Nagarscnker  (1972)  have  recommended  Z^  for  tests  of  hypotheses  I)  to 
III)  and  studied  the  general  distribution  problem.  The  density  of  has 
been  obtained  in  various  forms  by  Sudjana  (1973b)  and  further  by  Hart  and 
Money  (1976). 

Harmonic  mean  criteria  were  proposed  by  Pillai  in  1955  and  obtained 
approximations  to  the  densities  in  the  three  cases  (Pillai,  1955).  The 
exact  distribution  problem  and  the  moments  of  the  H?s^  criteria  have  been 
studied  by  Troskie  (1971,  Troskic  and  Money,  1972)  in  terms  of  zonal 
polynomials. 

Bagai's  statistic.  Y^,  has  been  studied  by  Bagai  in  several  papers 
(1962a,  1964b,  1965b,  1967,  1972a,  1972b)  giving  exact  and  limiting 
distributions  in  integral  forms  first  and  then  hypergeometric  functions  using 
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inverse  Mellin  transform.  Explicit  results  are  given  for  s*2(l)8  in  the 
exact  case  ami  s«5fl)10  in  the  Uniting  cast.  Further,  Consul  (1964b)  has 
obtained  the  limiting  distribution  for  5*20)8,  and  exact  distribution  in 
integral  form  using  inverse  Mellin  transform.  Further,  Mathai  (1970a, b,  1972b) 
has  obtained  the  exact  distribution  in  terns  of  Mcijcr's  G-function  through 
inverse  Mellin  transform.  Some  results  on  series  expansions  of  G-functions 
arc  also  given.  Further,  Oliver  (1972)  has  also  considered  the  distribution 
problem  in  terms  of  the  G-function. 

In  the  one-sample  case,  the  distribution  of  the  L statistic  for  the  test 
of  IV)  has  been  obtained  by  Nagarscnker  and  Filial  (1973)  and  percentage  points 
given  for  p*2(l)10.  Korin  (1968)  has  expressed  the  null  distribution  of 
-21og  L in  the  form  of  an  asymptotic  series  of  centra]  chi-squared  distributions 
and  computed  percentage  points  but  his  tables  arc  incomplete  for  small  v and 
p*3(l) 10  (Pearson  and  Hartley  (1972)).  Davis  (1971)  has  expressed  the 
percentage  points  of  -2  log  L in  terms  of  chi-squared  percentiles  using  a 
Comish-Fisher  inversion  of  Box's  series  but  his  tables  arc  also  incomplete 
in  regard  to  small  values  of  v for  the  values  of  p he  has  considered  i.e. 
p-b  and  10.  An  asymptotic  expansion  of  the  distribution  of  -2  log  L up  to 
order  v"3  has  been  given  by  Sugiura  (1969b)  involving  chi-souared  tenas 
inverting  characteristic  function. 

The  sphericity  criterion  of  Mauchly  (1940a, b)  for  test  of  V)  has  been 

considered  by  several  authors.  The  exact  distribution  of  W has  been  obtained 

by  Mauchly  (1940b)  for  p«2  and  by  Consul  (1967b)  for  p»3,4  and  6 and  further 

by  Consul  (1969),  Mathai  (1970a, b,  1971b)  and  Mathai  and  Rathic  (1970)  and  by 

John  (1972)  for  p»3.  The  expression  given  by  Consul  (1969)  is  in  tenas  of 

Meijer'sG-functions  (Mcijcr, 1946a,  1946b)  where  as  those  of  Mathai  and  Rathie 

arc  in  series  form.  Nagarscnker  and  Piliai  (1972,  1973)  have  obtained  the 

exact  distribution  of  W in  series  form  by  methods  developed  similar  to  the 

one  used  by  Box  (1949)  (also  sec  Andersen,  1958)  and  Nair  (1938,  1940)  and 

tables  of  5 and  1%  points  for  p*2(l)10  made.  Bagai  (1965b,  1972a)  has  also 

obtained  the  null  distributions  of  W as  a multiple  series.  A n asymptotic 

-2 

expansion  of  the  distribution  of  the  criterion  -2p  log  W to  order  v 
involving  chi-square  terms  has  been  given  by  .Anderson  (1958,  p.263),  where 
p is  a correction  factor  depending  upon  p and  v. 

Wilks  (1932)  defined  the  determinant  of  the  covariance  matrix  as  a 
scalar  measure  of  scatter  in  a multivariate  distribution  and  called  it 
generalized  variance.  Wilks  (1960)  has  discussed  the  relationship  of 


scatter  to  various  other  problems.  Hence  the  distribution  of  the  determinant 
of  the  SI1  n.ttrix  of  .1  sample  from  a p-variate  normal  population  has  been 
studied  !»>  varum*  author..  Ilu*  ratio  of  the  detemmants  of  the  -SI'  matrix  .uni 
the  covariance  matrix  can  he  expressed  as  a product  ol  p independent 

chi-squared  variables.  Wijsman  (1957)  has  studied  the  distribution  of  sample  ! 

’I 

generalized  variance.  Bagai  (1962b,  1965a)  has  obtained  the  distribution  in 
j the  null  case  (actually  in  the  noncentral  linear  case)  for  p«2(l)10  expressing 

| in  terms  of  integrals.  Consul  (1964a)  has  also  obtained  the  exact  distribution 

I * in  the  null  (linear)  case  using  inverse  Mellin  transform  and  special  cases  given 

for  p»2(l)7.  Again  Mathai  (1970a)  and  Mathai  and  Ratnic  (1971a)  have  obtained 
the  exact  distribution  in  the  null  (linear)  case  using  inverse  Mellin  transform 
and  G-functlons  and  expansions  by  calculus  of  residues,  psi  and  zeta  functions. 

Mathai  (1970a)  has  given  the  distribution  of  the  product  of  p independent 
gamma  variables  as  a G-function  and  similarly  that  of  the  product  of  p 
independent  beta  variables.  He  deduces  several  interesting  examples  from 
these  results  including  the  generalized  variance  for  the  gamma  case,  Wilks' 
criteria  and  others  for  the  beta  case.  Steyn  (1967)  and  Oliver  (1972)  have 
also  considered  the  distribution  problem,  the  latter  in  terms  of  inverse 
Mellin  transform. 

Khatri  (1967),  and  Pillai  and  Al-Ani  (Pillai,  Al-Ani  and  Jeuris,  1969) 
have  considered  a hypothesis  $ Ej  » E2,  6 > 0 unknown  which  includes  I)  as  a 
special  case  and  suggested  the  use  of  ratios  of  the  ch.  roots,  b./b^  (or 
gj/gj  in  the  limiting  case)  for  tests  in  this  connection.  The  latter  authors 
have  also  obtained  the  distribution  of  the  ratio  of  g^/gj  for  s»2,3  and  4. 

Krishnalah  and  Waikar  (1971a,  1972)  have  also  obtained  the  distribution  of 
the  ratios  of  successive  roots  and  each  root  to  the  maximum  root  in 
connection  with  simultaneous  tests  for  equality  of  latent  roots  against 
certain  alternatives  in  the  one-sample,  two-sample,  MANOVA  and  canonical 
correlation  cases.  Further,  Krishnaiah  and  Schuurmann  (1974)  have  obtained 
the  distribution  of  the  ratios  of  the  individual  roots  to  the  trace  in  the 
one-sample  case  and  some  percentage  points  obtained  in  connection  with  certain 
simultaneous  tests  of  the  hypothesis  (Krishnaiah  and  Waikar,  1971a,  1972). 
t In  this  connection  use  has  been  made  of  a relationship  established  by  Davis 

(1972c)  between  the  Laplace  transforms  of  the  ratios  of  individual  roots  to 
[ the  sum  of  the  roots  and  the  densities  of  the  above  individual  roots  in  the 

^ one-sample  case. 
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Further,  elementary  symmetric  functions  (esf)  of  the  eh.  roots  have  been 

of  some  interest  as  statistics  for  the  different  tests  of  hypotheses  of  which 

Hotelling's  and  Pi  Hal's  traces  arc  special  cases.  (Sec  Lemma  S.l)  Monotonicity 

properties  of  these  criteria  in  connection  with  test  1)  have  been  discussed  by 

Anderson  nnd  Oas  Gupta  (1964b) . Mijarcs  (1961)  has  shown  some  properties  of 

completely  homogeneous  symmetric  functions  and  certain  determinantal  results 

to  give  an  inverse  derivation  of  the  moments  of  V^.  The  method  has  been 

further  extended  to  the  moments  in  general  of  esf  of  the  b.  illustrating 

r r\  1 

explicitly  for  the  second  moment  of  V;  J . Mijarcs  (1964b)has  obtained  deter- 

minantal  expressions  for  moments  illustrating  the  results  to  third  moment 

of  second  esf  and  product  moments  of  first,  second  and  third  esf's.  Pillai 

(1965b),  in  addition  to  Lemma  5.1,  has  shown  some  relations  between  the  r**1 

moment  of  ufs|  „ „ and  ufs£  „ as  well  as  the  rth  moment  of  vfs]  _ and  a 
s-i  ,m,n  i gn-Tim+r 

r*\ 

linear  function  of  the  moments  of  order  up  to  r of  U,  J and  obtained 

l ,n,m*r 

approximate  upper  5 and  1%  points  for  n for  various  values  of  m an  n. 


Such  approximate  percentage  points  have  been  obtained  by  Dotson  (1968)  for 

vo5«  In  the  one  sample  ease,  the  second  esf  has  been  studied  by  Pillai 
2,ra,n 

and  Gupta  (1967)  in  regard  to  first  four  moments,  suggesting  an  approximation 
to  its  distribution. 


PART  II 

NOM-NIII.I.  DISTRIBUTION'S 


.s.  1NTROIHKTION 


This  part  of  the  paper  deals  with  the  non-central  distributions  of  the  eh. 
roots  in  connection  with  the  three  hypotheses  but  as  special  cases  of  Pillai's 
distribution  of  the  eh.  roots  under  violations  (Pillai,  197S).  The  exact,  sub- 
asymptotic  and  asymptotic  distributions  of  the  different  test  statistics  are 
considered  and  the  one-sample  case  also  treated.  Section  9 reviews  the  distri- 
bution of  the  ch.  roots  under  violations  covering  the  three  standard  distributions 
and  the  one-sample  case  and  in  addition,  treats  the  multivariate  beta  distribution. 
Further,  sub-asymptotic  distributions  arc  discussed  briefly  in  Section  10  and 
the  distribution  of  individual  ch.  roots  in  Section  11  dealing  with  both  exact 
and  asymptotic  cases.  The  distribution  problem  of  Hotel ling's  trace  is  considered 
next  in  Section  12  discussing  the  moments,  the  exact  distribution  and  asymptotic 
expansions.  Section  13  surveys  similar  studies  for  Pillai's  trace  and  Section 
14  of  Wilks'  criterion.  As  in  the  null  case,  other  statistics  are  treated  in 
Section  15,  namely,  Wilks'  statistic,  Z^,  Pillai's  harmonic  mean  criteria, 

Bagai's  statistic,  Y^,  and  the  It  test  of  equality  of  two  covariance  matrices. 
Actually,  these  critoria  except  the  harmonic  means  fall  under  a class  proposed 
by  Pillai  and  Nagarsenker  (1972).  Finally,  a short  critical  review  is  attempted 
in  Section  16  which  shows  that  in  spite  of  all  the  efforts  of  hundreds  of 
researchers  in  the  field,  the  results  achieved  so  far  generally  remain  un- 
satisfactory for  various  reasons. 
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•>.  non-mh.i.  ni.vruiiurnoKS  m Till  cn.  roots 

Hu*  jmnt  d i st rl luit ion  uT  the  eh.  roots  h»v  the  same  form  (see  hq  (2.1)  or 
(2.0)  or  (2.7)),  under  the  three  null  hypotheses  1)  to  III).  However,  in  the 
non-null  case,  the  joint  distributions  are  different  in  the  three  cases.  James 
(1964)  has  listed  them  systematically  and  the  reader  may  be  referred  to  his 
paper  for  a close  study.  Nevertheless,  in  the  next  section,  by  a suitable 
approach,  a single  expression  for  the  non-null  distribution  is  given,  of 
which  the  noncentral  distributions  for  I)  and  II)  fall  as  special  cases  and 
III)  also  follows.  In  this  section,  a few  preliminaries  will  be  discussed. 

It  has  been  shown  (Roy,  1957)  that  the  joint  density  of  the  ch.  roots  of 
S.S^1  for  I)  involves  as  parameters  only  the  ch.  roots  of  E.E^1  which  will  be 
denoted  by  A ■ diag(Xj,,...,Xp).  Such  result  is  true  for  II)  and  III)  as  well 
with  the  respective  parameter  matrices  ji  E'^yw'  having  ch.  root  matrix  (say) 

-1-1  - - - 2 9 9 

u « diagfwj wp)  and  having  ch.  root  matrix  P ■ diag(p“. . ,p“) . 

Here  j,  E^yii'  is  the  nonccntrality  matrix  in  the  nonccntral  Wishart  distribution 

of  S for  II)  (Anderson,  1946,  Anderson  and  Girshick,  1944,  Wcibull  } 1953,  James, 
*• 

1955a,  Her 2,  1955,  Constantine,  1963),  and  for  III),  the  covariance  matrix 
P ^ii  **i2 

E « (“  “ ) . (Further,  refer  Roy,  1957,  Anderson,  1958,  to  see  how  the  test 

q ^2  -22 
P 

of  II)  and  that  of  general  linear  hypothesis  lead  to  the  same  distribution 
problem) . 

The  method  employed  for  a unified  approach  to  the  distribution  problem  is 
to  make  the  parameter  matrix  partially  random  (denoted  "random"  hereafter) 
which  implies  diagonalization  by  an  orthogonal  transformation  II  and  integration 
over  H;  in  other  words  putting  a llaar  prior  on  H (see  James,  1964,  Constantine, 
1963,  Anderson,  1958,  for  Haar  measure)  leaving  the  latent  roots  of  the 
parameter  matrix  non-random.  Now  a theorem  is  given  below  (Pillai,  1975) 
concerning  the  distribution  of  the  characteristic  roots  of 

when  Sj(pxp)  has  a non-central  Wishart  distribution  with  df  and 
noncentrality  matrix  f.  E”*yy'  and  covariance  matrix  Ej,  and  SjCpxp)  has  an 
independently  distributed  central  Wishart  distribution  with  Vj  df  and  covariance 
matrix  E where  > p.  Let  F = diag.  (fp  . . . , where  fj,...,fp  are  the 

ch.  roots  of 


2* 


Theorem  0.1.  The  joint  density  of  is  given  by 

ts.U  c,(plvl,v,)e"wri?  |a|  s 1 1f|*}i  ♦ i f\  y 2 (fi“V 

- 1 - x>  j 

£ IrtCv^v,)^  Ck(XF(1  ♦ A F)_1)/k! 
k«0  k 

\ l\  4 C5(-»-‘a-‘)>.“  («)/»*, ^(1)1.0  « f,  if,  i ...  i f, 

d*0  6 * 

where  A is  "random",  A is  a positive  real  number,  m^^-p-l)  and 
c2(p,vrv2)  -In^p“ rp C v x * v2d ] / 1 rp i ) rp (f v2) rp * 

• •<  ^ k !>/«  1/!  im  f f 5*  fK#»  ?nnnl  nolvnomial  of 


where  Fp(a) 


ri  r(a-A(i-l)).  C (A)  is  the  tonal  polynomial  of 
i«l  ' * - 


degree  k corresponding  to  the  partition  « • (kj,...,kp)  of  k into  not  more  than 

p parts,  kj  ik,>  0,  k{  kp  - k.  The  generalized  Laguerre 

polynomial  lJ(A)  is  defined  in  Eq  (14)  of  Constantine  (1966)  and  a^fi  are 

constants  (scc~Eq  (20)  of  Constantine  (1966),  given  in  (9.3)  below, 

tabulations  in  Pillai  and  Jouris  (1969)).  Further,  (a)  ■ » (a-^(i-l))^ 

i«l 

with  (a)^»a(a+l) . . . (a+k-1) . 

It  may  be  noted  that  Cjfo.v^Vj)  - C(p,m,n)  given  in  the  null  case  and 
r (a)  is  proportioned  to  the  value  of  the  Wishart  integral  having  2a  df  with 
£ * I.  Zonal  polynomial  C^(A)  is  a homogeneous  symmetric  polynomial  in  the 
eh.  roots  of  the  symmetric  matrix  A such  that  (tr  A)»  E^C^A)  and  the  value  of 
C (I)  is  given  in  Eq  (21)  of  James  (1964).  James  (1M1*  1964)  has  defined 
Z (A)  by  the  relation  C^CA)  ■ [xjzK] ^5)  [2^k!/2k!)Z^(A)  where  X[2<]^  *S 
dimension  of  the  representation  [2k]  of  the  symmetric  group  on  2k  symbols 
(see  Eq  (19)  of  James  (1964)).  James  and  Parkhurst  have  tabulated  coefficients 
of  the  zonal  polynomials  expressed  both  in  terms  of  elementary  symmetric  functions 
of  the  ch.  roots  and  in  terms  of  sums  of  powers  for  k«l(l)12  and  are  given  in 
Harter  and  Owen  (1974).  Z J$)  « k!ak,  where  aJ(  is  the  k esf  in  the  ch. 
roots  of  A and  orthogonality  relations  to  obtain  other  zonal  polynomials  of 
degree  k are  given  in  James  (1964).  James  (1968)  has  also  given  an  alternate 
method  of  construction  of  the  zonal  polynomials  by  use  of  the  Laplace- 
Beltrami  operator.  Crowther  and  Young  (1974)  have  suggested  a third  method 
of  construction.  Further, 


i 


19. 2)  L^CAVHf.Vj)*  CKCI)}  - XC-l}ta<  TCT(A)/[(fv1jT  CT(I)]. 

The  a . coefficients  have  been  studied  in  addition  to  Constantine  (1966) 

* » i 

and  Pillai  and  Jouris  (1969),  by  Binghnm  (1974),  Muirhead  (1974)  and  others. 

The  density  of  bj,...,bp  can  be  obtained  from  (9.1)  by  the  transformation 
bj  ■ Af1/(1+Af^),l*l,.,.,p.  formula  (9.1)  yields  the  following  special  cases: 

(a)  For  ft«0,  from  Eq  (21)  of  Constantine  (1966),  (see  also  Hert  (1955)  page  487), 
l."(0)  * (fv.)  C (I).  Substituting  L?(0)  into  (9.1)  and  making  use  of 
(Constantine,  1966)  ^ 

(9.3)  C ( 1 -A) /C  (I)  - l $(-1)*  a C (A)/C  (I) , 

* - - * - ^aQ  j K I 1 ' - 

the  non-null  distribution  of  Khatri  (1967)  for  I)  is  obtained. 

b)  By  letting  A«I  and  A*1  and  by  using  the  relation  (9.2)  with  t«6»k,  the 

non-central  distribution  of  Constantine  (1963)  for  II)  is  deduced. 

Further,  if  ft  in  (9.1)  is  considered  a completely  random  matrix 
f Ej^'  M Y Y*  M'  eT^'  where  M(pxq)  is  a parameter  matrix  and  Y Y1  has  a central 
Wishart  distribution  N(q,v.,E.),  then  we  get  the  following  theorem:  (Pillai, 

1975) 

Theorem  9.2.  The  joint  density  of  f.,...,f  , where  ft  is  completely  random,  is 

ip  - 


given  by 


-*v.  -fv_  _ -^(v.+v,) 

(9.4)  C(p,v,v)|A|  1 1 1 + ft,|  ‘ 3|Fp|l*AF|  H (f  -f  ) 

£ 1 ^ - - - 1 - - - X J 

■ k 

I l lrt(Vj*v2))  C (XF(I*XF)"l)/k!j  l l a C (-X“1A“1)/C  (I) 
(>0  k 1 *■  * K - - - d«0  6 ,T  - 0 ~ 


k«0  k 
d 


tI0  \ w\x 


where  ft 


1 


i .1  x 

Ej  M'E^M  Ej. 


0 < f j £ f 2 ~ * * • = fp  < » * 


-1, 


In  (9.4)  taking  (I  ♦ ftp  ft^  = Pw,  A « I,  * vi+v2  an<*  the 

non-null  distributions  of  Constantine  (1963)  for  III)  is  obtained. 

The  non-null  distributions  of  Khatri  (1967)  for  I)  and  Constantine  for 
II)  and  III)  have  been  used  in  the  two-roots  case  to  compute  exact  powers  of 
the  criteria  1)  to  4)  by  Pillai  and  Jayachandran  (1967,  1968),  of  5)  by 
Sudjana  (1973)  and  of  individual  ch.  roots  by  Pillai  and  Dotson  (1969)  and 
Pillai  and  Al-Ani  (1970).  The  expressions  in  (9.1)  have  been  used  by  Pillai 
and  Sudjana  (1975)  to  study  exact  robustness  of  1)  to  4)  against  non-normality 
for  I)  and  against  the  violation  of  the  assumption  of  equal  covariance  matrices 


in  II).  Further,  Pillai  and  Hsu  (1975)  have  used  (9.4)  to  study  the  exact 

robustness  of  1)  to  4)  against  non-normality  for  III).  Zonal  polynomials  of 

degree  k ^ 0 have  been  used  in  the  above  computations. 

In  the  one-sample  case,  the  joint  density  of  the  eh.  roots  of  a Wishart 

matrix  in  the  non-central  case  with  unknown  E and  ft  has  not  been  worked  out. 

« *» 

But  the  density  when  a)  't  a o and  h)  E is  known  can  be  obtained  as  special 

m * 

cases  of  the  following  one  obtained  _rom  the  joint  density  of  B ■ diag. 

_i  - 

(b.,..,b  ) ■ AF(I*XF)  by  putting  <1  * v,B  and  making  v tend  to  infinity: 

Theorem  9.5.  The  joint  density  of  gj,..,gp  as  a limiting  ease  from  (9.1)  is 
given  by  (writing  v for  and  E for  \) 


(9.5J  K.  (p, v)etr*2 1 ^ | ^|5f^(v"I>"1)e“tr5  n (g.-g. 

i>j  1 * 


l l (CJGJ/k!)  I I a,  ACs(-X’1r1)Lf(v"P"1)(u)/(^v)4C“(I), 


k-0  < 


d-0  6 


0 < Si  m 4 8n  < "• 

±2  1 P 

where  Kj  (p,v)  » /[rp($v)rp(J,p)] . 

Note  that  Kj(p,v)  ■ K(p,m)  expressed  in  different  notations.  Formula  (9.5) 
yields  the  following  special  cases: 

a)  ft  ■ 0.  As  in  (9.1)  when  ft  ■ 0 applying  (9.3),  the  form  of  the  density 
given  by  James  (1960)  and  Pillai  and  Al-Ani  (Pillai  and  Al-Ani,  1967,  Pillai, 
Al-Ani  and  Jouris,  1969)  is  obtained. 

b)  E known.  Again,  as  in  (9.1),  putting  E ■ I,  X«l,  using  the  relation  (9.2) 
and  t«6«k,  the  form  given  by  James  (1961,  1964)  is  deduced. 

In  the  derivation  of  the  densities  above  for  the  eh.  roots,  there  are 

two  results  which  are  very  important  and  are  given  below. 

Lemma  9.1.  Let  R(pxp)  be  a positive  definite  symmetric  (pds)  matrix  and 

H(pxp)  he  an  orthogonal  matrix  with  positive  first  column  such  that 

H'R  II  » diag(r, , . . ,r  ) , where  r,,..,r  arc  the  eh.  roots  of  R.  Under  this 

---  1 P 1 P 

transformation  the  volume  element  dR  becomes  (James,  1954) 

• P 

(9.6)  dR  « n (r.-r.)  n dr.  dll 
- i>j  1 3 i-1  1 ” 

where  the  measure  dH  is  that  derived  by  the  exterior  product  of  differential 

forms  on  the  orthogonal  group.  With  this  measure  / dH  * 2pn*'p  /f  (£p), 

0(p)  “ P 

where  0(p)  denotes  the  group  of  all  orthogonal  (pxp)  matrices,  (James,  1954). 


3) 


Lemma  9.2  Let  S be  a symmetric  matrix  and  R,  a pds  matrix.  Then  (James, 
1960) 

(9.7)  / C Ill'S  II  R) (dll)  * t*  (S)C  (KJ/C  (U, 

0{p)  - - - * * * 

where  (dll)  is  the  invariant  llaar  measure  on  the  orthogonal  group,  normalised 

m 

so  that  the  measure  of  the  whole  group  is  unity. 

In  the  special  case  (b)  of  (9.1),  the  use  of  the  above  two  lemmas  will 
give  the  densities  of  F ■ Sl^'S.Sl^  and  B • (S,*S  ) "^S.  (S.+S,)“^  as  below. 


Theorem  9.4  The  density  of  F is  given  by  (James,  1964) 

(9.8)  C1(p,v1,v2)itrli|F|a|UF|  ‘V1+V2)iFi^(v,-v2)^vi;«F(UF)-1) 


and  the  density  of  B is  given  by 

(9.9)  C^p.v^vpe’^lBf  |l-B|n1F1(f.(v1+v2);^,nB), 

where  n«f.(v,-p-l)  ,Cj  (p,  Vj , v2)  « rpr^(v1+v2)]/(rp(^v1)rp(§v2)) 

and  hypcrgeometric  function  of  matrix  argument  (James,  1964,  Constantine,  1963) 

P*VC  i » • • »cp»di  * * • *dq*^)  " k]0  I (d^Jl.Td^  T T“* 

The  density  (9.8)  is  known  as  the  non- central  multivariate  F distribution 
and  (9.9)  non-centraimultivariate  Beta  distribution.  Kabe  (1963)  derived  (9.9) 
when  rank  of  fl  is  2 and  Sitgreaves  (1952)  for  the  two  roots  case. 

If  Vj  <?  p £ v0  instead  of  p <.  v^,v0  as  assumed  earlier,  the  densities  in  this 
case  can  be  derived  from  those  above  by  making  the  following  substitut&cns: 

(9.10)  (v^v,,?)  (p.v^v^-p.Vj). 


10.  SUB-ASYMPTOTIC  DISTRIBUTIONS  OF  THE  CH.  ROOTS 

In  recent  studies  on  sub-asymptotic  distributions  of  ch.  roots,  the 
approach  has  generally  been  to  maximize  an  integral  of  the  form 

I«J  _P_(ci  * • * • *c_»di » * • • »dn»  AM  R,,,)(dll) 

0(p)P  ^ 1 P 1 

where  0 (p)  is  the  group  of  orthogonal  matrices  H(pxp)  with  respect  to  which 

the  maximization  is  carried  out,  A * diag(a. , . . . ,a  ) , R * diag(r. ,...,r ), 

*•  * P - * P 

(dH)  is  the  invariant  or  Haar  measure  over  the  group  0 (p)  normalized  so  that 

the  measure  of  the  whole  group  is  unity,  F^  is  a hypcrgeometric  function  of 

matrix  argument  and  cj» • • »Cp»d^» « « • »dq  are  functions  of  df  and  are  positive 

real  numbers.  In  the  one-sample  or  covariance  matrix  case,  G.  A.  Anderson  (1965) 


has  shown  that  the  maximum  of  the  integrand  with  p»q»0,  R,  the  sample  eh. 

root  matrix  and  A”*,  the  matrix  of  population  eh.  roots  assumed  to  be  distinct, 

% 

is  attained  when  II  « diag.  (+1 , . . ,-t  1) . James  (1969)  has  extended  the  study  to 
the  ease  uf  an  extreme  multiple  ]>opulation  root.  Chang  (1970)  and  Li  and 
Pillai  (Li  and  Pillai,  1970,  Li,  Pillai  and  Chang,  1970)  have  found  the  form 
for  ti  the  same  as  in  the  one-sampic  case  when  maximizing  the  1*0  1‘  tsgrand  in 
the  two-sample  (two  covariance  matrices)  problem  with  p»l  and  q»0.  Chattopadhyay 
and  Pillai  (1970,  Chattopadhyay,  Pillai  and  Li,  1976)  have  generalized  the 
results  of  the  previous  authors  not  only  to  pF^  hypergcometric  function  but 
also  when  a^'s  are  equal  in  each  of  several  sets. 

The  asymptotic  expansion  for  large  df  is  obtained  by  evaluating  the 
integral  for  small  neighborhoods  of  the  matrices  I!  ■ (±1,...,±1).  While 
the  earlier  authors  considered  the  asymptotic  expansions  of  a single 
covariance  matrix  or  two  covariance  matrices,  Chattopadhyay  and  Pillai 
(1971b,  1973,  Chattopadhyay,  Pillai  and  Li,  1976)  have  extended  the  study  to 
MANOVA,  jFj,  and  canonical  correlation,  ,Fj,  cases  and  from  one  extreme 
multiple  root  situation  to  that  of  one  multiple  root,  extreme  or  intermediate. 

In  extending  the  work  further  to  the  case  of  several  multiple  population  roots, 
the  method  used  by  James  (1969)  was  not  found  to  be  suitable  in  view  of  the 
fact  that  the  invariance  of  a function  with  respect  to  the  choice  of  a submatrix 
in  the  orthogonal  matrix  used  there,  docs  not  extend  to  the  similtaneous 
invariance  with  respect  to  choices  of  several  submatrices  as  is  needed  to 
extend  the  method.  Chattopadhyay  and  Pillai  (1971,  1973,  Chattopadhyay,  pillai 
and  Li,  1976)  have  shown  that  the  integrand,  if  A ■ diag.(a1,...,ar,ar+1,*!,' 

%•*!’  ’ * * ’*r+d‘  *^,ar+d  ^ is  0PtiBUW  if  " dia*,CI0W»  t?l » * * * *?d^  * where 
In(r)  « diag. (±1, . . . ,±1) , H.(t.xt.),i»l,...,d,  arc  orthogonal  matrices.  It  is 

maximum  or  minimum  depending  on  the  ordering  of  a/s. 

The  method  of  expansion  is  to  consider  the  transformation  H ■ exp[S]  where 

S(pxp)  is  a skew  symmetric  matrix  of  the  form 

/ c (0)\ 


S-l  S(i)  I S(0)(rxp);S(i)  * (s{i},  , S^),  Sj(i)  (t.xfr+tj*.  .n.^)) , 

- Vd>/  ~ 

$2  (t^xt^)  * 0,  S j (t^x (t^+|+ . . +t j) ) , i«l , . . . ,d-l ; 

and  S(d)  « (s[d) , S2(d)),  s{d3(tdx(p-td)),  S2(d)(tdxtd)  « 0. 


r 


With  the  appropriate  ordering  of  a^*s,  the  integrand  is  maximized  when  H 

is  of  the  above  form  and  for  large  df  v,,  the  whole  integral  is  concentrated 

•» 

around  its  unique  maximum,  lienee  there  is  no  loss  of  generality  in  using  the 

above  transformation  provided  the  constant  factor  is  adjusted  as  at  least  one 

maximizing  set  is  covered  by  this  substitution.  (See  Pillai  and  Chattopadhyay, 

1970,  1973).  For  large  v,  and  r^'s  and  a^'s  well  spaced,  most  of  the  integrand 

will  be  given  by  small  values  of  S and  terms  of  the  expansion  can  be  computed 

* 

in  view  of  the  transformation  II  ■ cxp[S]. 


Constantine  and  Muirhead  (1976}  have  very  recently  studied  the  asymptotic 

expansion  for  the  distribution  of  the  eh.  roots  in  the  one-sample  ease  for 

unknown  t and  ft«0,  two-sample  F-case  for  I)  and  type  I case  fo t II)  for  eh. 

•»  « * 

roots  of  ft  Urge.  Their  method  takes  the  above  sub-asymptotic  approach  as 
* 

starting  point  to  derive  the  first  term  of  the  expansions  but  then  use  a 

partial  differential  equation  approach  (Muirhead  1970a  , Constantine  and 

Muirhead,  1972)  to  successively  calculate  the  other  terms  up  to0(v’3), 

-3  -2 

0(v,*v.)  and  Q(v  ♦v,)  respectively. 

In  the  one-sample  case,  Muirhead  and  Chikuse  (197S)  have  obtained  an 

asymptotic  expansion  in  terms  of  normal  density,  of  x^  « (v/2)^(t^-l) 

i>l,..,p,  using  the  sub-asymptotic  expansion  of  Anderson  (1965),  where 

aj's  are  the  ch.  roots  of  E“*.  The  terms  up  to  order  v”1  have  been  obtained. 

Bingham  (1972)  has  employed  a paramatrization  of  the  rotation  group 

0 (p)  of  pxp  orthogonal  matrices  with  determinant  *1  in  terms  of  their  skew 

symmetric  parts  to  derive  for  p*3  an  explicit  expansion  for  qFq(Z,H) , a 

hypergeometric  function  of  two  matrix  arguments  appearing  in  the  distribution 

of  the  ch.  roots  of  a pxp  tfishart  matrix.  On  the  basis  of  a numerically 

derived  simplification  of  the  low  order  terms  of  this  series,  an  asymptotic 

expansion  for  p«3  of  qFq  in  terms  of  products  of  ordinary  confluent 

hypergeometric  series  is  conjectured.  Limited  numerical  exploration 

-8 

through  terms  of  order  v has  indicated  that  the  new  series  is  several  orders 
of  magnitude  more  accurate  than  the  scries  from  which  it  was  derived. 

11.  INDIVIDUAL  CH.  ROOTS 

The  study  of  individual  ch.  roots  in  the  non-central  case  has  been 
carried  out  by  several  authors.  Zonal  polynomial  series  have  been  useful 
here. Pillai  (1966c, 1970)  has  obtained  the  exact  cdf  of  the  largest  root 
for  II)  and  III)  for  the  two-roots  case  and  three-roots  linear  case  and 
computed  powers  of  the  test.  Pillai  and  Jayachandran  (1967)  haveobtained 
an  approximation  to  the  largest  root  in  two-  } three-  and  four-roots  linear 


I 
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cases.  Again,  Pillai  and  Jayachandran  (1968)  ha*  obtained  the  cdf  in  the 
two-roots  case  and  powers  computed.  Pillai  and  botson  (1969)  have  obtained 
the  cdf  of  the  individual  ch.  roots  for  II)  and  III)  for  two  and  three-roots 
cases  and  Pillai  and  Al-Ani  (1970)  for  1)  ind  power  computed  in  all  cases. 

In  all  the  above  study  the  tonal  polynomials  up  to  sixth  degree  have  been  used. 
Al-Ani  (1970)  has  obtained  the  distributions  of  the  second  largest  roots  for  1) 

II)  and  III)  in  terms  of  zonal  polynomials.  Further,  Pillai  and  Al-Ani  (1969) 
have  obtained  the  distributions  of  the  smallest  and  second  smallest  roots  for 

I) ,  II)  and  III)  and  the  one-sample  case.  In  the  MANOVA  case  i.e.  hypothesis 

II) ,  Hayakawa  (1967)  and  independently  Khatri  and  Pillai  (1968a)  have  obtained 
the  density  of  the  largest  root  in  a beta  function  series  form.  The  density 
involves  two  types  of  coefficients,  gt  _ and  b seen  from  the  following 

1C  | T K jT 

relations: 

(11.1)  CJA)Ct(A)  « £ g^.C6(A), 

where  k ■ (kj,..,kp),  t ■ (tj,...,tp)  and  6 ■ (dj,...,dp)  such  that  d ■ k*t, 

where  < is  a partition  k k,  x of  t and  6 of  d.  Further,  if  M.  ■ diag(M,l), 

» * «• 

where  M ■ diag(m1,...,mp_l),mi's  being  the  ch.  roots  of  M,  then  (Hj)  could 
be  written  as  ^ 

(11.2)  C (M  ) « l l bK  C (M) . 

-*  t-0  T K*T  T * 

Khatri  and  Pillai  (1968a)  have  tabulated  g*  coefficients  for  d < 7 and 

k,t  . 

b coefficients  for  k < 6.  Hayakawa  (1967)  has  tabulated  h,  coefficients 

£ Kp T 

which  are  g times  a factor, for  d < 4.  Similarly  Hayakawa  has  a coefficients 
corresponding  to  Khatri  and  Pillai's  b 's.  Further  Khatri  (1967)  has  obtained 
the  density  of  the  largest  root  for  I)  as  a jF,  hypergeometric  function. 

ai  and  Sugiyama  (1969)  have  obtained  simpler  power  series  expressions 
fc  .ue  density  and  the  cdf  of  the  largest  root  for  the  MANOVA  case  than 
obtained  before  and  also  derived  those  of  the  largest  root  for  test  of  I)  and 

III)  both  in  the  type  I and  type  II  cases.  De  Waal  (1969a)  has  also  considered 
the  distribution  of  the  largest  root  for  III).  Further^Troskie  and  Money  (1972) 
have  obtained  the  distribution  of  the  smallest  root  for  III).  Waikar  (1973)  also 
obtains  the  joint  density  of  f^  and  fp  as  well  as  gj  and  g^.  Khatri  (1972)  has 
obtained  the  distribution  of  the  largest  as  well  as  the  smallest  roots  for  I),  II) 
and  III)  in  exact  finite  series  form  involving  zonal  polynomials  when  n^t^-p-l) 
is  a non-negative  integer.  Based  on  Khatri' s finite  series,  Venables  (1973)  has 
presented  an  algorithm  for  the  numerical  evaluation  of  the  null  distribution  of 
bp.  Krishnaiah  and  Chattopadhyay  (1975)  have  discussed  the  extensions  of  the 


i 
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i 
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pfaffian  method  to  the  non-central  case. 

In  the  one-sample  case,  Hayakawa  (1.969)  has  obtained  the  density  of  the 

largest  root  of  the  non-central  Wishart  matrix  with  known  t as  a series  in 

* 

Hcrmite  polynomials.  For  unknown  I but  Q>0,  Khatri  (1972)  has  given  the 
distribution  of  the  largest  root  as  well  as  the  smallest  root  in  finite 
series  of  zonal  polynomials  for  non-negative  integral  values  of  m. 

The  non-central  distributions  of  Khatri  for  I),  II)  or  III)  or 
one-sample  case  although  in  finite  series  form  is  of  limited  use  since  they 
are  generally  extremely  complex.  Hence  Pillai  and  Saweris  (1974b)  have 
obtained  new  forms  based  on  Khatri's  results,  as  finite  or  infinite  series, 
but  useful  for  further  work  and  generally  more  rapidly  convergent  than  those 
obtained  earlier  by  Pillai  and  Sugiyama  (1969).  They  have  also  given  some 
approximations  to  these  distributions. 

In  regard  to  asymptotic  distributions,  Hsu  (1941a)  has  obtained  for  II) 
the  following  theorem: 

Theorem  11.1.  Let  zt  ■ N*(2wj>  4wh)  *(fi-wh),(i»mh_1+l,.^;h«l,..,t),zi«Nfi, 


(i>r*l,...,s).  Then  the  limiting  distribution  of  Zj,..,z  as  the  sample  sizes 
tend  to  infinity  is  given  by  the  density  (sample  ch.  roots  ordered  in  descending 
order  of  magnitude) 

* fay* J > • • • » 


f (*  J # • • » * 


) . • . f(z  , , • . . 

"i  Vi  1 


where 


' u . u y x^ 

f(x, x ) - 2"*u(  n r(*i))_1  n (x.-xje  1 , 

i-l  i<j  1 3 


> X1  m"  »xu  > 
•.  s-r 


flCzr>l'“'rs)  * 2‘*(p‘r)(l’1"r)n*(s’r)[  n r(iq-^r-fi^)r(ii))-1 

S 1,1 

n n (z.-z.)  n z icq-*-1)*-*  f xi  , • > z £...£  z > o, 
i«r+l  j»i*l  1 3 i«r*l  1 imT +1  T *l  5 

where  w^  >...>  wt  > 0 are  the  ch.  roots  of  fi  of  multiplicity 

respectively  and  mQ«0,  m^jj,  ♦..♦  j^(h«l,.. ,t) ,r«mt  and  q>max(p,.l-l) 

and  N is  the  total  of  t sample  sizes. 

Hsu  (1941b)  has  obtained  similar  asymptotic  distributions  for  III)  namely 
canonical  correlation.  Further,  Anderson  (1948)  has  adopted  Hsu's  proofs  to 
obtain  the  limiting  distributions  of  the  ch.  roots  in  the  regression  problem 
'(linear  hypothesis)  and  also  that  of  the  ch.  roots  of  a sample  covariance 


j 
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matrix  in  the  non-central  case.  In  all  cases,  it  is  interesting  to  note  that 

for  different  h the  ch.  roots  are  asymptotically  independently  distributed. 

Anderson  huu  further  studied  (1951,  1963)  in  the  one-sample  case  the 

limiting  distribution  of  the  ch.  roots  and  vectors  in  view  of  principal 

component  analysis.  If  the  characteristic  roots  of  l are  different,  the 

■* 

deviations  of  the  ch.  roots  of  the  sample  covariance  matrix  from  the 
corresponding  population  quantities  are  asymptotically  independently 
normally  distributed  with  variance  2X?/v;i»l ,. . . ,p.  Further,  when  some  of 

the  ch.  roots  of  51  arc  equal,  the  asymptotic  distribution  of  sample  ch.  roots 

* 

and  vectors  also  have  been  studied. 

_2 

Muirhead  (1970b)  obtained  an  asymptotic  expansion  up  to  order  v for  the 

largest  root  in  the  one-sample  case  based  on  his  expansion  for  jFj  function. 

But  it  is  of  limited  interest  because  it  is  valid  only  over  the  range  of  values 

less  than  the  smallest  eh. root  of  E. Further,  MuiThead  and  Chigusc  (1975)  using 

* 

the  system  of  pde's  of  Muirhead  (1970a)  have  derived  an  asymptotic  expansion 

up  to  order  v"1  of  the  distribution  of  x » (v/2)^(t  a -1)  and  siml  rly  that 

P P P 

of  Xj.  They  have  also  obtained  the  marginal  density  of  xl  from  the  asymptotic 
joint  distribution  discussed  in  the  previous  section  in  terms  of  normal  densities 
which  agrees  with  the  result  of  Suglura  (1973). 

12.  HOTELLING'S  TRACE 


The  non-central  distribution  of  Hotelling's  trace  has  been  obtained  by 
Constantine  (1966)  through  inverse  Laplace  transform  in  the  following  form: 

(12.1)  C3(p,v1#v2)(U(pJ)  1 l C(-UCp}  /(ipv^^k.'H^CVj+v,))^^), 


k-0 


where  C3(p,vJ,v2)  « C2(p,Vj  .VjjTp^v^/r^pVj) , 

m « ^(v^-p-i) ,L^(H)  is  the  generalized  Lagurrc  polynomial  given  in  (9.2)  and 

p < Vj,v,.  The  series  is  convergent  for  [U^l  <1.  The  density  of  U(p^  for 

V1  < P i v2  can  be  obtained  from  (2.1)  using  (9.10).  The  central  case  can  be 

obtained  by  putting  f(«0  as  in  special  case  (a)  of  (9.1). 

% 

Pillai  (1973,  Pillai  and  Sudjana,  1974)  has  used  (12.1)  to  suggest  the 
following  form  for  the  density  of  U^: 

(12.2)  C,(p,v.,v,)(Ul',))iPVl’I/(l*UtP)/p)4P<Vl>V2)  I «-l)kBk(Ut!’)/p)l7(l*U(P,/p)1 

v ' .k-1  k-1  k'° 

where  U.  = tpk/(f.pv.)vk!]I(f.(v  ))  L™(n)-  l [ Jl  (Kvj+v^p+jj/fk-r)  1)E ,Ej*0. 

K 1 * * 1 *■  K K ' r=o  j=r 


For  p»2  and  fi«0,  (12.2)  reduces  to  the  exact  form  given  by  Hotelling 
(Constantine,  1966).  P'or  higher  values  of  p no  verification  has  been 
possible. 

PI 

The  exact  distribution  of  Uv~'  using  up  to  the  sixth  degree  zonal  ^ 
polynomial  has  been  obtained  by  Pillai  and  Jayachandran  (1967,  1968)  and 
used  for  computing  powers  of  the  appropriate  tests.  Khatri  (1967)  has 
obtained  the  distribution  of  Hotelling's  trace  for  I)  following  Constc;tine 
(1966)  in  series  form  again  convergent  for  |U^|  < l.  Pillai  and  Sudjana 
(196S)  have  further  extended  the  results  of  Constantine  and  Khatri  into  a 
single  density  function  starting  from  (9.1),  Pillai  and  Hsu  (1975)  for  canonical 
correlation  using  (9.4)  und  specialized  for  the  two-roots  case  following  Pillai 
and  Jayachandran,  for  robustness  studies  against  non-normality  for  1)  and  III) 
and  against  the  violation  of  the  assumption  of  equality  of  covariance  matrices 
for  II). 

In  regard  to  the  study  of  moments  for  II),  Hsu  (1940a)  h«s  obtained  the 

first  two  moments  of  although  the  second  moment  is  slightly  in  error  as 

shown  by  Khatri  and  Pillai  (1967)  who  have  obtained  the  first  four  moments. 

Earlier,  Ghosh  (1965)  has  obtained  the  first  moment  of  and  the  variance 

for  values  of  p-3  and  4.  Pillai  (1964c,  1966c)and  Khatri  and  Pillai  (1965,  1966) 

have  obtained  the  first  four  moments  of  in  the  linear  case  based  on 

independent  beta  variables  which  extend  to  the  linear  case  and  an  approximation 

to  its  distribution  suggested  based  on  the  first  three  moments.  Further,  from 

the  study  of  the  first  four  moments,  Khatri  and  Pillai  (1968b)  have  given  two 

approximations  to  the  distribution  of  in  the  general  non-cectral  case. 

Constantine  (1966)  has  derived  the  general  moment  of  in  the  form 

(12.3)  E(U(p))k  ■ (-l)k  r L"(-Sj)/(*(p*l-v2))r,v2  > 2k+p-l. 

< 

In  regard  to  asymptotic  distributions  of  considerable  work  has  been 
carried  out  by  various  authors.  Asymptotic  expansion  of  the  distribution  of 
Hotelling's  T?  * v,U^  has  been  obtained  by  Ito(1956)  in  the  null  case  for  II) 

-2 

up  to  order  v2  . The  non-null  distribution  was  given  by  Siotani  (1957)  and 
later  by  Ito  *.1960)  up  to  order  v^1.  Further,  Siotani  (1968,  1971)  extended 

_ o 

his  result  up  to  order  v,“.  The  above  authors  generally  used  Taylor  expansion 

and  perturbation  techniques  of  James  (1954)  either  on  the  cdf  or  characteristic 

functions.  Fujikoshi  (1970)  has  given  two  methods  of  obtaining  the  non-null 
■>  _2 

distribution  of  Tq  up  to  order  v?  by  using  hypergeometric  function  and  Lsguerre 
polynomial  of  matrix  argument  using  characteristic  function  and  Laplace  transform. 


Both  lend  to  the  same  result  which  is  the  same  as  that  of  Siotani  (1968).  The 
above  expansions  all  are  generally  in  terms  of  non-central  ehi -squared  variables. 
Further,  Muirhead  (1972b)  has  also  derived  an  expansion  up  to  order  which 
can  be  readily  util iced  to  give  further  terms  in  the  expansion  and  which  may  be 
difficult  to  obtain  by  the  use  of  earlier  methods,  ilayakawa  (1972)  has  also 
given  an  asymptotic  expansion  of  the  distribution  of  T"  to  order  v,”  by 

W k 

obtaining  some  formulae  for  generalized  Laguerre  polynomials  of  matrix 
argument  and  univariate  Laguerre  polynomials. 

For  III),  Fujikoshi  (1970)  has  obtained  an  asymptotic  expansion  for  Tq 
for  testing  the  hypothesis  of  independence  between  two  sets  of  variates  under 
sequence  of  alternatives  converging  to  the  null  distribution  with  a rate  of 
convergence  v“y(y  ^0).  By  utilizing  a close  relationship  between  II)  and  III), 

* “ o 

the  asymptotic  expansion  of  the  distribution  of  Tq  up  to  order  has  been 
studied  for  Y"l.  Fujikoshi  has  also  computed  the  powers  using  the  asymptotic 
distribution  for  II)  for  the  three-roots  case.  Further,  Lee  (1971a)  has 
extended  Ito's  result  for  II)  to  III).  The  asymptotic  expansions  in  all  the 
work  above  arc  generally  involving  non-central  chi-squared  terms. 

Again,  Sugiura  and  Nngao  (1973)  have  obtained  asymptotic  formulae  for  the 
distribution  of  TQ  for  III)  up  to  order  (vi*v2)  *n  terms  of  normal  distribution 
function  and  its  derivatives.  Method  is  based  on  differential  operators  on 
symmetric  matrix  by  Siotani  (1957),  Ito(1960)  and  others. 

For  I),  if  S.  and  S,  are  defined  as  in  (9.1)  except  that  fi«0, 

Chattopadhyay  and  Pillai  (1971a)  have  obtained  asymptotic  expansions  for  cdf 

i -1  _1 

and  percentile  of  T2  to  order  v«  . The  expansions  hold  when  E.El  * I+E  and 

II  ^ * ek  - - 

|ch.  Ej  <1,  i«l,..,p,  where  ch. (A)  denotes  the  i ch.  root  of  A.  Taylor 
expansion  and  perturbation  techniques  have  been  employed.  Pillai  and  Saweris 
(1973)  have  extended  the  results  of  Chattopadhyay  and  Pillai  to  order  and 

also  included  terms  involving  f-^f^  in  v”*  order  terms  which  were  neglected  by 

the  earlier  authors  where  is  the  (i,j)th  element  of  EjE^1.  Further,  Pillai 
and  Saweris  (1974a)  have  extended  the  work  of  Chattopadhyay  (1972)  who  derived 
an  asymptotic  expansion  up  to  order  v2  for  the  cdf  and  percentile  of  T‘  with 
Sj  and  S2  as  in  (9.1)  with  fl|0.  The  extension  was  to  include  f.^  f^  terms 
neglected  by  Chattopadhyay.  Hsu  and  Pillai  (1975a)  have  extended  these 

results  further  to  include  terms  up  to  order  but  neglecting  terms  of  higher 
powers  in  some  Taylor  expansions  to  keep  the  number  of  terms  to  a manageable 
level  for  computation.  In  all  this  work  perturbation  technique  is  used.  The 
asymptotic  expansion  obtained  for  the  non-centrul  case  has  been  extended  also  to 


the  canonical  correlation  case  by  making  ii  random.  Numerical  results  of  powers 
have  been  computed  from  the  expansions  by  Pillai  and  Sawcris  (1974)  and  Hsu  and 
Pillai  (1975)  for  robustness  studies  for  1)  and  III)  against  non-normality  and 

11)  against  the  violation  of  assumption  of  equality  of  covariance  matrices. 

13.  PILLAI 'S  TRACE 

The  exact  distribution  of  is  more  difficult  to  obtain  in  the 
non-ccntralcasc  than  even  that  of  II ^ for  various  obvious  reasons  and  it 
has  not  been  derived  so  far.  The  exact  distribution  for  the  two-roots  case 
has  been  obtained  by  Pillai  and  Jayachandran  (1967,  1968)  for  I)  II)  and  III) 
using  up  to  sixth  degree  tonal  polynomial.  Khatri  and  Pillai  (1968a)  have 
obtained  the  density  of  for  II)  as  tonal  polynomial  scries  but  is 
convergent  only  for  jv^|  < 1. 

In  regard  to  the  moments,  Pillai  (1964c,  1966c)and  Khatri  and  Pillai 

(1965,  1968a)  have  obtained  the  moments  of  in  the  linear  case  based  on 

the  idea  of  independent  beta  variables.  Again, Khatri  and  Pillai  (1967)  have 

shown  that  the  moments  of  can  be  obtained  from  moments  of  for  II) 

where  r < p is  the  rank  of  fl.  The  first  two  moments  of  have  been  worked 

out.  Further,  Pillai  (1968)  has  obtained  the  mgf  of  for  I)  II)  and  III) 

but  the  mgf  of  I)  has  been  obtained  by  Khatri  (1967)  with  a slight  error. 

The  Tcsults  are  in  zonal  polynomial  series  form  which  involve  the  a „ 

",  * 

coefficients  in  (9.3)  which  have  been  tabulated  only  for  k $ 8 (Constantine, 
1966,  Pillai  and  Jouris,  1969). 

Pillai  and  Sudjana  (1975)  have  extended  the  work  of  Pillai  and 

Jayachandran  (1967,  1968)  for  robustness  studies  for  I)  and  II)  using  (9.1) 

and  Pillai  and  Hsu  (1975)  for  III)  using  (9.4)  as  in  the  case  of  U^. 

They  Have  also  extended  (Pillai  and  Sudjana,  1972)  the  mgf  of  Pillai  (1968) 

starting  from  the  density  (9.1)  and  Pillai  and  Hsu  (1975)  starting  from  (9.4). 

In  regard  to  asymptotic  distributions,  Fujikoshi  (1970)  has  obtained  an 

asymptotic  expansion  of  chi-squared  terms  for  the  distribution  of  for  II) 
_2 

up  to  order  v2  , by  using  Pillai 's  mgf  (1968)  and  using  certain  formulae  for 
weighted  sums  of  zonal  polynomials.  He  has  shown  the  usefulness  of  the 
asymptotic  expansion  by  obtaining  approximate  upper  5 and  1%  points  and 
comparing  with  those  of  Pillai  and  Jayachandran  (1967)  and  by  computing  powers 
of  the  test  for  three-roots  case.  Similarly  for  III),  Fujikoshi  (1970) 
has  obtained  an  asymptotic  expansion  using  Pillai* s mgf.  For  y*l  (see  Section 

12)  , the  expansion  consists  of  non-central  chi-squared  terms  and  has  been 

_2 

derived  up  to  order  . Further,  Lee  (1971b)  has  obtained  an  asymptotic 
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expansion  for  to  order  and  powers  computed  and  upper  percentage 

points  of  V<P>  also  computed  by  Hi  1 1 -Davis  method  (1968).  Lee  (1971a) 

also  has  extended  to  III)  the  expansion  for  II).  Again  Sugiura  and  Nagao 

(1972)  have  obtained  asymptotic  formulae  for  the  distribution  of  for 

III)  up  to  order  in  terms  of  normal  distribution  function  and  its 

derivatives  using  the  method  of  differential  operators  on  symmetric  matrix 

by  Siotani  (1957),  Ito  (1960)  and  others.  Fujikoshi  (1972)  has  further 

-2 

obtained  an  asymptotic  expansion  to  order  Vj  for  II)  by  employing  a method 
similar  to  that  of  Sugiura  and  Nagao(1971).  The  asymptotic  formulae  are 
given  in  terms  of  chi-squared  variables. 

14.  WILKS'  CRITERION 

The  exact  non-central  distribution  of  Wilks'  criterion  has  been  derived 

by  Pillai,  Al-Ani  and  Jouris  (1969)  for  I),  II)  and  III)  using  inverse  Mellin 

transform  in  terms  of  Mcijer's  G-functions.  Pillai  and  Nagarsenker  (1972) 

P . k 

has  obtained  the  distributions  of  a statistic  n b. (1-b.)  of  which 

i-1  x 1 

Pillai-Al-Ani-Jouris  results  are  special  cases.  Pillai  and  Sudjana  (1975) 

have  extended  the  results  starting  from  (9.1)  and  Pillai  and  Hsu  (1975)  using 

(9.4).  Earlier,  Pillai  and  Jayachandran  (1967,  1968)  have  obtained  the  exact 

distribution  of  Wilks'  criterion  for  I),  II)  and  III)  in  the  two-roots  case 

and  computed  powers  in  each  case.  Das  Gupta  (1972)  has  shown  that  the  null 

and  the  non-null  linear  distribution  of  Wilks'  criterion  have  the  monotone 

It  property.  Das  Gupta  and  Perlman  (1973)  have  shown  that  the  power  of  W^ 

strictly  decreases  with  the  dimension  and  Vj  in  the  linear  case. 

Further,  in  view  of  independent  beta  variables  apynri  H'  which  extends  to 

the  linear  case  with  the  non-centrality  parameter  inv  1*. ,-d  r J.y  in  the  density 
2 

of  tjj  in  (6.2)  as  can  be  seen  from  (9.9)  (or(9.8))  b>  <ran-> formation  as  in 
Section  6,  Gupta  (1971)  obtained  the  density  of  Wilks'  ri.rtion  for  II)  in  the 
linear  case  for  number  of  roots  2 to  5 using  the  convolution  technique  (see 
Section  6).  Asoh  and  Okamoto  (1969)  have  expressed  the  criterion  as  a product 
of  conditional  betas  in  the  non-null  case  while  Kabe  and  Gupta  (1972)  have 
given  an  explicit  procedure  for  decomposing  the  non-central  multivariate  beta 
density  in  terms  of  densities  of  independent  beta  variables  extending  the  work 
of  Pillai  (1964c,  1966c)and  Khatri  and  Pillai  (1965).  Again,  Hart  and  Money 
(1975) "have  obtained  the  distribution  of  Wilks'  criterion  in  the  linear  case 
in  a general  form  and  have  developed  an  algorithm  for  calculating  exact  powers 
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and  percentiles.  Also,  tabulations  of  powers  anti  percentiles  have  been 
carried  out  for  selected  number  of  variates  and  df.  Trettcr  and  Walster 
(1975)  also  have  obtained  the  distribution  in  the  linear  case  as  mixtures 
of  beta  distributions.  Nagarsenker  (1976)  has  derived  the  distribution  of 
I),  II)  and  III)  in  the  general  case  in  simpler  computable  forms  as  mixtures 
of  beta  distributions  as  well  as  mixtures  of  gamma  distributions  which  are 
more  useful  than  Pillai-Al-Ani-Jouris  forms  involving  G-functions.  The 
G-functions  are  not  easy  to  evaluate  in  all  cases.  Mathai  (197G'a)  has  given 
scries  expansions  for  certain  cases  of  G-functions  and  shown  the  use  of  the 
series  for  distributions  of  Wilks'  criterion  derived  by  Pillai-Al-Ani-Jouris. 
Mathai  (1973a)  has  brought  about  the  points  involved  in  the  evaluation  of  the 
G-function  in  terms  of  computable  expressions. 

In  regard  to  approximate  or  asymptotic  distributions,  Ito(1962)  has  made 
numerical  power  comparisons  of  test  of  II)  based  on  T2  and  Wilks'  criterion  in 
large  samples  in  the  linear  case.  Mikhail  (1965)  has  given  such  comparison  by 
an  approximate  method  using  the  lower  order  moments  based  on  and  W^. 

further,  Posten  and  Bergmann  (1964)  have  given  asymptotic  expansions  for  II) 
using  characteristic  functions  in  the  linear  case  and  Roy  (1966)  has  also 
tabulated  the  power  throush  an  approximate  method.  Sugiura  and  Fujikoshi 
(1969)  have  derived  for  II)  an  asymptotic  expansion  involving  chi-squared 
terms  to  order  f°r  Wilks'  criterion  using  the  characteristic  function 
expressed  in  terms  of  hypergeometric  function  of  matrix  argument.  Olkin  and 
Si°tani  (1964)  have  shown  for  III)  that  the  limiting  distribution  of 
(n ')  W<P>  is  normal.  For  III),  the  h*^1  moment  of  the  criterion  has  been 
worked  out  first  and  then  studying  the  characteristic  function  Sugiura  and 
Fujikoshi  (1969)  have  also  obtained  an  asymptotic  expansion  up  to  order 
*n  t-crms  of  normal  distribution  function  and  its  derivatives. 

Fujikoshi  (1970)  has  computed  approximate  powers  based  on  these  expressions 
in  the  three-roots  case.  Lee  (1971a)  has  given  an  asymptotic  expansion  for 
III)  using  Sugiura-Fujikosbi  (1969)  expansion  for  II)  as  also  Sugiura  (1969). 
Power  comparisons  have  been  made  (Lee  1971b)  of  U^p\  and  W^  for  II) 
in  the  three-and  four-roots  cases.  Muirhead  (1972a)  has  obtained  an 
asymptotic  expansion  for  the  distribution  of  tre  up  to  order  (^♦v^-3  for 
local  alternatives  in  terms  of  non-central  chi-squared  distribution  and 
density  by  obtaining  a solution  of  the  system  of  pde's  satisfied  by  the  mgf. 

This  extends  the  work  of  Sugiura  (1969).  Pillai  and  Nagarsenker  (1972)  has. 
obtained  an  asymptotic  expansion  for  I)  up  to  order  (v^+Vj)  ^ in  terms  of 
normal  distribution  function  and  its  derivatives  and  Subrahmaniam  (1975)  has 
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extended  the  same  to  order  (Vj*v,)"*^“  using  pdo's  of  Muirhcad  (1970a)  giving 
power  tabulations  of  the  test,  the  actual  tr  test  and  test. 

15.  OTHER  STATISTICS 

As  in  the  null  case,  Wilks1  statistic,  Z^,  will  be  considered  first 

in  this  section.  The  approach  to  the  distribution  problem  is  generally  the 

same  as  that  for  Wilks'  criterion  but  there  has  been  less  interest  about  this 

study  in  view  of  the  inferior  nature  of  the  tests  regarding  power  compared  to 

the  more  outstanding  test  criteria  like  U^,  and  W^.  De  Waal  (1968) 

has  derived  the  non-central  density  of  Z^  in  terms  of  an  integral  expression. 

Pillai  andpNagarsenker  (1972)  have  obtained  the  distribution  of  a statistic  of 

the  form  II  b?(l-b,)^  of  which  Z^  is  a special  case.  The  distribution  has 
i-1  1 1 

been  derived  in  terms  of  H-functions  (Mathai,  1970a)  as  a result  of  inverse 
Mcllin  transform.  The  exact  distribution  has  been  obtained  for  I),  II)  and 
III)  and  asymptotic  expansions  have  been  obtained  to  order  (vj'fv2)”1  *n  tcrws 
of  normal  distribution  and  derivatives  for  I)  under  local  alternatives. 

.t  12 

Subrahmaniara  (1975)  has  extended  the  results  to  order  * Sudi*na 

(1973b)  has  obtained  the  general  form  of  the  distribution  of  Z^  starting 
from  (9.1)  and  (9.4)  and  obtained  the  distributions  for  I),  II)  and  III)  as 
special  cases.  He  has  also  given  power  tabulations  based  on  detailed  study 
of  the  two-roots  case  and  compared  the  powers  with  those  of  U^,  V^,  W^ 
and  largest  root.  The  distribution  in  the  non-central  linear  case  has  been 
studied  by  Troskie  and  Mcney  (1974)  following  Gupta  (1971)  obtaining  the 
distribution  of  Z^  for  II)  for  p«2,3,4  and  5.  They  have  also  extended  the 
results  to  III).  Further,  Hart  and  Money  (1976)  have  applied  their  (1975) 
algorithm  for  Wilks'  criterion  with  modification  to  obtain  non-central  linear 
percentiles  and  powers  of  Z^  and  carried  out  numerical  power  studies  whose 
findings  agree  with  those  of  the  extensive  study  of  Sudjana  (1973b).  Dc  Waal 
(1968)  has  obtained  an  asymptotic  distribution  for  Z^  assuming  that  n is  a 
fixed  matrix.  Fujikoshi  (1972)  has  derived  an  asymptotic  expansion  up  to 
order  involving  chi-squared  terms  by  inverting  the  characteristic  function 
expressed  in  terms  of  hypergeometric  function  with  matrix  argument. 

Regarding  the  three  Harmonic  mean  criteria,  l/P,\»l,2,3,  following 
Constantine  (1966)  and  Khatri  and  Pillai  (1968a),  the  densities  have  been 
studied  by  Troskie  (1971)  which  are  convergent  only  in  the  range  |Hp^  | <1  . 

The  densities  have  been  obtained  in  the  context  of  III).  Further,  Troskie  and 
Money  (1972)  have  obtained  the  densities  for  II)  and  mgf's  for  II)  and  III). 
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The  exact  distribution  of  Bagai’s  statistic,  , has  been  derived  by 

Pi  1 Ini  and  Nagarsenker  (1972)  in  terns  of  G-  and  li-functions  using  inverse 

Mcllin  transform  for  I),  II)  and  III).  Sugiura  (1969a)  has  obtained  an 

usymptotic  expansion  for  I)  for  large  sample  sizes  and  local  alternatives. 

In  regard  to  the  lr  test  of  I),  it  falls  as  a special  case  gf  the 

criterion  of  Pillai  and  Nagarsenker  (1972)  given  above,  namely,  It  bT (1 -b . )** 

i-1  1 1 

with  a « Vj/2  and  b * v^/2.  The  exact  distribution  has  been  given  by  the 

authors  in  terms  of  G-  and  H- functions.  Further  they  have  obtained  asymptotic 

expansions  for  the  criterion  inverting  the  characteristic  function  in  two 

_2 

forms,  the  first  as  chi-squared  series  up  to  order  *n<*  *econ<*  as 

chi-squared  series  up  to  order  ra  “ where  ra»p(v^v2)  where  p is  a correction 
factor  (Anderson,  1958,  p.255)  but  both  under  local  alternatives. Subrahmaniam 
(1975)  have  extended  these  expansions  using  the  pde's  of  Muirhead  (1970a)  to 
order  (vi+v2)”3  an<*  m”3  respectively. 

In  the  one-sample  case,  the  exact  non-central  distribution  of  the  L 
statistic  has  not  received  a great  deal  of  attention.  Sugiura  (1969b)  has 
obtained  asymptotic  expansion  of  the  distribution  of  the  criterion  -21og  L 
up  to  order  v-1  involving  normal  distribution  function  by  inversion  of  the 
characteristic  function. 

In  regard  to  the  sphericity  criterion,  Girshick  (1941)  has  obtained  the 
non-null  distribution  of  sphericity  criterion  in  the  two-roota  case  in  series 
form.  Pillai  and  Nagarsenker  (1971)  have  obtained  the  exact  distribution  in 
terms  of  zonal  polynomials  and  G- function  by  inverse  Mel 1 in  transform.  Khatri 
and  Srivastava  (1971)  have  also  derived  the  distribution  in  similar  form. 
Steffens  (1974)  has  obtained  power  tabulations  for  the  two-roots  case  obtaining 
an  F series  form  for  the  distribution  equivalent  to  that  of  Girshick 
(Girshick,  1941,  James,  1966)  and  shows  that  the  power  of  the  sphericity  test 
excels  that  of  a bivariate  t test  of  sphericity  suggested  in  the  paper  except 
in  the  close  neighborhood  of  the  hypothesis.  Muirhead  (1976)  has  also  computed 
powers  of  W in  the  two-roots  case  based  on  the  non-null  distribution  which  he 
has  obtained  as  mixtures  of  F distributions  similar  to  that  of  Steffens. 

As  regards  asymptotic  distributions,  Gleser  (1966)  has  given  a normal 
asymptotic  form  and  so  has  Sugiura  (1969b).  Sugiura  (1969b)  has  given  the 
asymptotic  expansion  of  -2  plogW  to  order  v"1  involving  normal  distribution 
function  and  its  derivatives  by  inverting  the  characteristic  function. 
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Sugiura  (1972)  has  also  obtained  the  Uniting  distribution  of  a locally  best 
invariant  test  discussed  by  John  (1972)  and  Sugiura  (1972)  given  by 
tr  G“/(tr  G)2. 

As  regards  generalized  variance,  Mathai  (1970a,  1972a)  has  obtained  the 

cxuct  non-central  distribution  of  the  dcteminant  of  the  SP  matrix  as  series 

involving  zonal  polynomials  and  G-functions  obtained  through  inverse  Hell in 

transform.  Further,  the  G-functions  are  evaluated  in  computable  series  form. 

Fujikoshi  (1968)  has  obtained  the  Uniting  distribution  when  and  has 

•»  % 

derived  an  asymptotic  expansion  of  the  distribution  of  the  sample  generalized 

-3/2 

variance  to  order  v when  fi  is  constant.  Further,  Fujikoshi  (1970)  has 
extended  the  study  to  the  case  n « b obtaining  the  expansion  up  to  order 

_ t/') 

v in  terms  of  normal  distribution  functions  and  its  derivatives  by 
inverting  the  characteristic  function.  Sugiura  and  Nagao  (1971)  also  have 
obtained  a similar  expansion. 

In  regard  to  elementary  symmetric  functions,  Khatri  and  Pillai  (1968b) 

have  obtained  a recurrence  relation  to  find  the  moments  of  the  i1^  esf  in  p 

bj's  as  well  as  p f^'s  in  the  non-central  linear  case  in  terns  of  those  of 

partitioned  submatrices.  Pillai  and  Jouris  (1969)  have  given  exact  expressions 

for  the  first  two  moments  of  _ and  first  three  moments  of  uP^,  „ based  on 

i,»tn  i|i»n 

the  results  of  Khatri  and  Pillai  (1968b)  and  also  the  third  moments  of  the 

second  esf  of  a matrix  in  the  non-central  means  case  (James,  1961).  The  first 

two  moments  in  question  have  been  obtained  earlier  by  Pillai  and  Gupta  (1968) 

in  terms  of  generalized  Lagucrre  polynomials.  The  expected  values  of  the  esf's 

in  the  p ch.  roots  of  a Wishart  matrix  in  the  central  case  has  been  obtained 

by  De  Waal  (1972a)  and  a conjecture  in  the  non-central  case  which  was  proved 

to  be  correct  by  Shah  and  Khatri  (1974).  The  ht^'  moment  of  the  trace  of  a 
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non-central  Wishart  matrix  with  Z«c  I has  been  studied  by  Nel  (1971)  and 

first  six  moments  explicitly  computed  for  general  Z (1972).  De  Waal  (1973) 

and  De  Waal  and  Nel  (1973)  have  further  studied  the  expected  values  of  the  esf's 

in  the  Wishart  and  correlation  matrices  cases.  Further,  Fujikoshi  (1970)  has 

obtained  an  asymptotic  expansion  of  the  distribution  of  the  trace  of  a 

non-central  Wishart  matrix  under  the  assumptions  D«0(1)  and  0(v)  to  order 
-3/2 

v in  terms  of  normal  distribution  function  and  its  derivatives  by 
inverting  the  characteristic  function. 


16.  A SHORT  CRITICAL  REVIEW 


Over  the  past  fifty  years,  considerable  work  has  been  promoted  in 
multivariate  distribution  theory  and  a large  number  of  papers  written  in  the 
area  has  been  devoted  to  characteristic  root  distributions.  While  the  null 
distribution  of  the  ch.  roots  was  obtained  as  early  as  1939,  the  non-null  ch. 
root  distributions  were  derived  only  in  the  sixties  and  those  too  in  series 
forms  which  are  not  immensely  useful  in  view  of  convergence  difficulties  and 
the  inability  to  compute  the  general  terms.  However,  even  more  unsatisfactory 
has  been  the  study  of  the  exact  distributions  of  the  test  criteria  which  are 
functions  of  the  ch.  roots  and  which  have  been  established  to  have  various 
optimum  properties  (see  Section  3).  The  exact  distributions  of  individual 
roots  have  been  obtained  explicitly  for  small  number  of  roots  but  when  the 
number  of  roots  becomes  larger,  one  has  to  resort  to  approximations.  Such 
approximations  have  been  suggested  only  for  the  largest  (smallest)  roots 
(see  Pillai,  1954,  19S6a,  1965,  Section  4).  Attempts  to  obtain  approximations 
for  other  roots  have  failed  (see  Davis,  1972a,  Section  4).  The  study  of  the 
exact  null  distributiontof  the  trace  statistics  has  met  with  less  success. 

Nanda  (1950),  Pillai  and  Jayachandran  (1970),  Davis  (1970b),  Krishnaiah  and 
Chang  (1972),  and  others  have  worked  on  the  exact  null  distribution  of 
obtaining  limited  results  for  small  number  of  roots  and  small  values  of 
arguments  m and  n.  However,  what  they  have  been  in  full  agreement  is  to 
recommend  the  use  of  Pillai' s tables  (1960)  of  approximate  percentage  points 
to  practitioners.  While  Pillai  has  given  such  tabulations  for  number  of 
roots  from  2 to  8,  Mijares  (1964)  has  extended  the  tables  up  to  number  of 
roots  50.  Again,  the  exact  null  distribution  of  Hotelling's  trace  has  been 
investigated  by  several  authors  with  somewhat  disappointing  results. 

Constantine  (1966),  Pillai  and  Young  (1971) , Pillai  and  Sudjana  (1974),  Davis 
(1968,  1970c)  and  others  have  worked  on  the  problem  but  the  general  non-null 
distribution  has  not  been  obtained.  Here  again,  the  practitioner  has  to  resort 
to  Pillai's  tables  (1960)  giving  approximate  percentage  points  for  number  of 
roots  2 to  8.  However,  progress  has  been  made  in  regard  to  the  exact  null 
distribution  of  Wilks'  criterion  in  view  of  inverse  Mellin  transform  approach 
which  was  introduced  by  Nair  (1938)  along  with  the  differential  equation 
approach  and  theory  of  residues,  which  are  some  of  the  important  approaches  in 
distribution  theory  followed  by  several  authors  (see  Section  6) . 
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Further,  approximations  to  the  distribution  of  Kilks'  criterion  were  available 
in  the  literature  for  a long  time  (Bartlett,  1938,  1947a,  Rao,  1948).  In 
addition,  the  exact  null  distribution  of  Z^  follows  that  of  W^,  and  those 
of  other  tr  test  criteria,  for  example,  L,  sphericity  and  others  also  have  been 
obtained  through  inverse  Mel  1 in  transforms.  The  study  of  the  distributions  of 
the  determinants,  for  example,  Bagai's  statistic,  generalized  variance,  has 
been  facilitated  by  the  use  of  G-functions  (see  Section  7). 

In  regard  to  the  non-null  distributions,  the  exact  distribution  problem 
has  generally  met  with  disappointing  results.  The  exact  distributions  of  the 
largest  (smallest)  or  other  roots  have  been  considered  by  many  authors 
(Pillai  and  Sugiyama,  1969,  Al-Ani,  1970,  Khatri,  1972,  Pillai  and  Saweris, 
1974b  and  others)  in  finite  or  infinite  series  forms  involving  zonal  polynomials 
whose  general  terms  cannot  be  explicitly  computed  and  with  slow  or  questionable 
convergence  in  the  infinite  case.  The  exact  non-null  distributions  of  or 

U(P>  have  not  been  obtained  in  the  general  case.  The  distributional  forms 
involving  zonal  polynomial  series  derived  in  this  connection  for  or 
arc  convergent  only  for  |V^|  or  |U^  | < 1 (Constantine,  1966 , Khatri  and 
Pillai,  1968a, Khatri,  1967  and  others).  However,  the  general  momehts  of 
has  been  obtained  (Constantine,  1966,  Khatri  and  Pillai,  1967)  and  the  mgf's 
of  (Pillai,  1968)  all  involving  zonal  polynomials  and  some  constants  whose 
general  expressions  are  unavailable.  The  exact  non-null  distributions  of  Wilks' 
criterion  have  been  derived  in  terms  of  G-functions  (Pillai -Al-Ani  and  Jouris, 
1969)  using  inverse  Mellin  transform.  The  G-functions  cannot  be  easily 
evaluated  in  all  cases.  But  computable  series  form  expressions  for  special 
cases  of  G-functions  involved  in  the  distribution  of  W^  and  others  have  been 
obtained  (Mathai,  1970a).  Attempts  have  also  been  made  to  give  the  distribution 
as  mixtures  of  beta  or  gamma  distributions  (Tretter  and  Walster,  1975, 
Nagarsenker,  1976).  As  in  the  null  case,  the  non-null  distribution  of  Z^ 
is  similar  to  that  of  W^  and  those  of  Bagai's  statistic,  sphericity  and 
generalized  variance  involve  G-functions  and  zonal  polynomials. 

While  a great  deal  remains  yet  to  be  accomplished  on  the  exact  null  and 
non-null  distributions,  extensive  work  has  been  carried  out  on  asymptotic 
distributions.  The  sub-asymptotic  expansions  treated  in  Section  10  are 
interesting,  however  the  distributional  forms  obtained  are  not  of  readily 
usable  form.  Starting  with  the  sub- asymptotic  expansions,  some  asymptotic 
expansions  have  been  obtained  in  usable  form  (Muirhead  and  Chiguse,  1975, 
Constantine  and  Muirhead,  1976,  Section  10)  but  the  sub-asymptotic  nature  has 
disappeared. 


Asymptotic  expansions  have  been  obtained  in  a large  nuaber  of  cases  by 
inverting  asymptotic  formulae  of  the  characteristic  functions  or  mgf  (I to, 
I960,  SJotani,  1971,  Sugiura  and  Fujikoshi,  1969,  Fujikoshi,  1970,  Huirhead, 
1970b  and  others).  In  one  set  of  expansions,  use  is  aide  of  soae  formulae 
for  weighted  sums  of  zonal  polynomials  and  generalized  laguerre  polynomials 
(Fujikoshi,  1970,  Sugiura,  1971)  in  another  set,  Taylor  expansions  and 
perturbation  techniques  are  used  (Ito,  1956,  Siotani,  1971,  Chattopadhyay  and 
Pillai,  1971,  Pillai  and  Saweris,  1973,  Hsu  and  Pillai,  1975  and  others)  and 
in  a third  set  de's  are  employed  (Davis,  1968,  1970a,  1970b,  1972a,  Huirhead, 
1970b  and  others)  to  obtain  solutions  for  mgf's  or  density  functions.  The 
results  have  been  obtained  generally  for  large  sample  sizes  in  terms  of 
chi-squared  distribution,  central  or  non-central,  or  normal  distribution. 

Some  expansions  in  the  non-central  case  are  for  local  alternatives. 

The  usefulness  of  these  asymptotic  results  have  been  very  little  so  far 
for  the  practitioner.  Some  results,  like  the  asymptotic  independence  and 
normality  of  the  ch.  roots  in  the  one-sample  case  for  principal  component 
analysis,  have  been  very  helpful.  But  then  the  unanswered  question  is: 

How  large  should  the  sample  size  be  to  use  these  results? 

In  conclusion,  it  has  to  be  recognized  that  quite  a large  number  of 
papers  have  been  written  in  the  field  of  ch.  roots  distributions  employing 
various  methods  of  approach.  Merits  and  demerits  of  some  of  these  methods 
have  been  discussed  by  Mathai  (1973b).  Although  considerable  efforts  have 
been  made  by  various  authors  on  exact  distribution  theory  especially  of 
different  test  criteria,  the  results  have  not  been  encouraging.  Again,  even 
though  the  exact  non-null  distributions  of  the  ch.  roots  have  been  available 
for  sometime,  it  seems  desirable  to  have  these  obtained  in  some  other  more 
useful  forms  than  the  existing  ones  i.e.,  series  involving  zonal  polynomials 
whose  convergence  is  slow  or  questionable  and  general  terms  are  not  com- 
putable. The  recent  progress  made  in  asymptotic  distribution  theory  remains 
yet  to  be  appreciated  and  assimilated  by  practitioners. 
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comparisons  made  In  connection  with  tests  of  three  multivariate  hypotheses. 
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covariance  matrix.  A brief  critical  review  has  also  been  attempted.  For 
convenience  in  organization,  the  review  has  been  given  in  two  parts:  Part  I. 


